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ABSTRACT 

The observed stellar mass function (SMF) is very different to the halo mass function 
predicted by ACDM, and it is widely accepted that this is due to energy feedback from 
supernovae and black holes. However, the strength and form of this feedback is not 
understood. In this paper, we use the phenomenological model GALFORM to explore 
how galaxy formation depends on the strength and halo mass dependence of feedback. 
We focus on "expulsion" models in which the wind mass loading, (3 is proportional 
to l/w^isk' 71 = 0, 1, 2 and contrast these models with the successful Bower et 

al. model (B8W7), for which oc l/v^ k . A crucial development is that our code 
explicitly accounts for the recapture of expelled gas as the system's halo mass (and 
thus gravitational potential) increases. While models with high wind speed and mass 
loading result in a poor match to the observed SMF, a model with slower wind speed 
matches the flat portion of the SMF between M* ~ 10 9 -10 n /j _1 M q . When combined 
with AGN feedback, the model provides a good description of the observed SMF above 
10 9 /i _1 M Q . In order to explore the impact of different feedback schemes further, we 
examine how the expulsion models compare with a further range of observational data, 
contrasting the results with the B8W7 model. In the expulsion models, the brightest 
galaxies are assembled more recently, and the specific star formation rates of galaxies 
decrease strongly with decreasing stellar mass. The expulsion models tend to have a 
cosmic star formation density that is dominated by lower mass galaxies at z = 1 — 3, 
and dominated by high mass galaxies at low redshift. These trends are in conflict with 
observational data, but the comparison highlights some deficiencies of the B8W7 model 
also. The experiments in this paper give us important physical insight to the impact 
of the feedback process on the formation histories of galaxies, but the strong mass 
dependence of feedback adopted in B8W7 still appears to provide the most promising 
description of the observed universe. 



1 INTRODUCTION 

A central issue in the study of galaxy formation is to under- 
stand the connection between the mass of galaxies and the 
mass of their associated dark matter haloes. This problem 
is not trivial: whilst the mass function of dark matter haloes 
predicted by the cold dark matter paradigm has a relatively 
steep slope (dlog N/dlogM ~ —0.9 over the range of haloes 
mass relevant to the galaxy formation), the observed mass 
function of galaxies is characterised by a shallow Schechter 
function. The abundance of galaxies is nearly independent 
of stellar mass over the interval M* ~ 10 8 - 10 10 ' 5 h' 1 M Q , 
whilst at greater masses the abundance of galaxies declines 
exponentially. This raises two fundamental questions: i) why 
is the stellar mass function so flat below ~ 10 10 ' 5 h~ 1 M@, 
when the abundance of haloes is such a strong function of 
mass, and ii) what physical processes are responsible for the 



exponential suppression of galaxies with mass greater than 
~ 10 10 ' 5 h _1 M (White & Frenk 1991; Benson et al. 2003)? 

In this paper we investigate the astrophysical processes 
that map the halo mass function to the galaxy stellar mass 
function (SMF). There are two common approaches to this 
problem. Arguably, the most appealing route is to numeri- 
cally integrate the set of differential equations that describe 
the rudimentary astrophysical processes (e.g. gravity, hydro- 
dynamics, radiative cooling, star and black hole formation) 
using a local framework (e.g. a set of particles or a grid). 
This approach aims to evolve the cosmological density fluc- 
tuation power spectrum reflected by the cosmic microwave 
background using an ab initio description of the problem. 
This system of equations commonly over-predicts the uni- 
versal abundance of stars — a problem commonly known 
as the 'over-cooling problem' (e.g. Katz, Weinberg, & Hern- 
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quist 1996; Balogh et al. 2001; Keres et al. 2009; Schaye et 
al. 2010). 

In order to get closer to a description ol the observed 
universe, an additional set of processes (collectively known 
as "feedback") must be introduced to couple the energy, 
mass and metals returned by supernovae (SNe) and black 
holes (AGN) to the surrounding gas. As well as accounting 
for the stellar mass and gas content of galaxies, feedback 
enriches the intergalactic medium with metals, potentially 
resulting in an orthogonal set of observational constraints. 
However, the wide range of scales involved in this prob- 
lem (from AU scales to tens of Mpc) makes it infeasible 
to model these processes from first principles, forcing re- 
course to phenomenological, or 'sub-grid', treatments (eg., 
Springel & Hernquist (2003); Springel, Di Matteo, & Hern- 
quist (2005); Okamoto, Gao, & Thcuns (2008); Schaye & 
Dalla Vecchia (2008); Booth & Schaye (2009)). 

Because the properties of model galaxies are remarkably 
sensitive to the details of sub-grid models, an alternative ap- 
proach is to establish a set of equations that describes the 
same astrophysical processes on macroscopic scales, typi- 
cally averaged over the physical scale of a galaxy. This ap- 
proach is adopted by phenomenological, or 'semi-analytic' 1 
models, such as the GALFORM code considered here (Cole 
et al. (2000); Bower et al. (2006), hereafter Bow06), see 
also Kauffmann, White, & Guiderdoni (1993); Hatton et 
al. (2003); De Lucia et al. (2006); Somerville et al. (2008); 
Guo et al. (2011)). For example, rather than computing the 
star formation rate at each point in a galaxy, such models 
typically compute the total rate of star formation over the 
entire galaxy, and assume that the rate is a simple function 
of global galaxy parameters (such as the gas mass, disk size 
and rotation speed) . This leads to an alternative set of dif- 
ferential equations, that provide a macroscopic description 
of the physics. A phenomenological model then solves these 
relatively simple equations within the merging hierarchy of 
structure formation. Macroscopic treatments are by defini- 
tion approximate, but their use is ubiquitous in all branches 
of physics. When applied carefully and within an understood 
range of validity, the resulting descriptions can lead to sig- 
nificant physical insight. 

In this paper, we apply the GALFORM model to seek 
an understanding of the mass-dependence of the equations 
of feedback from galaxies. We focus initially on the role 
supernova-driven winds play in establishing the galaxy stel- 
lar mass function on scales below ~ 10 10 ' 5 /i _1 M Q . Although 
strong suppression of galaxy formation in low mass haloes is 
clearly necessary, our understanding of the process remains 
incomplete, largely because of the difficultly of calculating 
the effects of feedback from fundamental physical princi- 
ples. The difficulty primarily arises because the interstellar 
medium in which SNe explode is inherently multiphase (and 
magnetized), and does not behave like an ideal gas. Thus, 



1 The term "semi-analytic" commonly used in the literature, but 
is missleading. With modern computing power it is no longer 
critical that the resulting equations can be solved analytically. 
The important distinction is that this type of model provides a 
macroscopic description of the relevant processes. This approach 
is key to gaining insight into the problem. Such models arc usually 
referred to as "phenomenological" in other areas of physics. 



accurately calculating the impact of even a solitary super- 
nova is extremely challenging, as the result depends strongly 
on the initial density of the medium into which energy is 
injected. As a sequence of SNe explode, a network of low 
density channels (or 'chimneys') is established that allow 
the supernova ejecta to escape into the halo of the galaxy 
(McKee & Ostriker 1977; Efstathiou 2000; de Avillez & Bre- 
itschwerdt 2007). However, without detailed calculations it 
impossible to estimate how the mass outflow rate and the 
specific energy of the outflowing material will depend on the 
star formation rate and the mass of the host galaxy. 

In view of this uncertainty, all currently feasible sim- 
ulations of galaxy formation parameterise the effect of su- 
pernovae and include it as a sub-grid calculation (Springel 
& Hernquist 2003; Dalla Vecchia & Schaye 2008). In hydro- 
dynamical simulations, winds are most commonly modeled 
by adding thermal energy to gas particles (or cells), or by 
giving gas particles a velocity kick. Although adding ther- 
mal energy seems a promising route, the energy injected is 
easily radiated away if the temperature of the particles is 
too low (Katz, Weinberg, & Hernquist 1996). Equally the 
injection of kinetic energy may drive thermodynamic shocks 
into the surrounding gas, and this thermal energy may also 
be radiated. These issues arise because the codes treat the 
complex multiphase interstellar medium (ISM) as a single 
fluid. Several schemes have been developed to circumvent 
the problem; one approach is to decouple the relevant par- 
ticles from the numerical scheme for a period of time, ei- 
ther by preventing cooling for a period (Brooks et al. 2007), 
or by decoupling the particles from hydrodynamical forces 
(Springel & Hernquist 2003; Okamoto, Nemmen, & Bower 
2008; Oppenheimer & Dave 2008). Particles therefore retain 
the energy injected by feedback for a period of time, eas- 
ing their escape from the dense interstellar medium of the 
galaxy. An alternative approach is to stochastically heat or 
kick the particles so that their energy is sufficiently high that 
their cooling time is long, or the shocks they drive are suffi- 
ciently strong that radiative losses are small (Dalla Vecchia 
& Schaye 2008; Creasey et al. 2011). 

Because of the difficulty in interpreting the multi-phase 
nature of the outflow, the mass loading and velocity of winds 
are not yet strongly constrained by observation (but see 
Martin (2005); Weiner et al. (2009); Chen et al. (2010); 
Rubin et al. (2011) for recent progress). Most hydrodynami- 
cal calculations have therefore adopted the simplest possible 
model, in which the mass loading and velocity of winds are 
independent of the system in which the feedback event is 
triggered. This also simplifies implementation of the feed- 
back, since there is no requirement to estimate the envi- 
ronment of the ISM on the fly, for example by determin- 
ing the mass of the dark matter halo, or the local gravi- 
tational potential. An exception is the momentum scaling 
model of Oppenheimer & Dave (2006), where the wind pa- 
rameters are set following a radiatively driven wind model 
(Murray, Quataert, & Thompson 2005), and a number of 
similar schemes implemented in the Overwhelmingly Large 
Simulations (OWLS Schaye et al. 2010). 

By contrast, most phenomenological ("semi-analytic") 
models of galaxy formation assume that the parameters de- 
scribing feedback adjust to ensure that the specific energy of 
outflows are matched to the binding energy of the halo. Con- 
servation of total energy thus ensures that the mass load- 
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ing of outflows is greater in dwarf galaxies than in larger 
systems (Dekel & Silk 1986). Such models are partially mo- 
tivated by arguments pertaining to the porosity of the in- 
terstellar medium: analytic models (e.g. Efstathiou 2000) 
consider the formation of channels in the multiphase ISM 
and suggest that the porosity of the ISM is self-regulating 
and determined by the gravitational potential of the disc. In 
general, phenomenological models further assume that ma- 
terial expelled from the disc is recaptured on a timescale 
proportional to the dynamical time. This coefficient is al- 
lowed to be significantly larger than unity in some models 
in order to approximate the effect of expulsion of gas from 
the halo. 

Current phenomenological models present a coherent 
picture for the formation of galaxies in a cosmological con- 
text, and provide an excellent explanation of many diverse 
data-sets. Although the GALFORM model has been de- 
veloped to explain the observed properties of galaxies (eg., 
Bow06, Font et al. (2008); Lagos et al. (2010)), it has also 
been shown to explain the X-ray scaling relations of groups 
and clusters (Bower, McCarthy, & Benson 2008, hereafter 
Bow08) and the optical and X-ray emission from AGN 
(Fanidakis et al. 2011). The models can be used to gener- 
ate convincing mock catalogues of the observable Universe 
(Cai et al. 2009) and applied to test the procedures used to 
derive physical parameters from astronomical observations, 
and to identify the priorities for the next generation of as- 
tronomical instruments. These successes have been driven 
by the inclusion of two key components of the models: i) 
galaxy winds that scale strongly with halo mass and ii) a 
"hot-halo" mode 2 of AGN feedback in which the cooling of 
gas in quasi-hydrostatic haloes is suppressed. By altering 
the parameterisation of the feedback schemes here, we in- 
vestigate whether the Bow06 choice is optimal, or whether 
alternative schemes can similarly reproduce the properties 
of observed galaxies. 

A major advantage of semi-analytic models is that they 
enable the effects of sub-grid parametrization to be explored 
quickly and easily (Bower et al. 2010). We exploit this aspect 
in this study, to explore the effect of various descriptions of 
outflows of material from galaxies. This requires us to gen- 
eralise the GALFORM model to include the possibility that 
gas is expelled from the potential of dark matter haloes. Pre- 
vious efforts to 'calibrate' the semi-analytic method against 
numerical simulations have included gravity, hydrodynamics 
and cooling, but not effective feedback (Benson et al. 2001; 
Helly et al. 2003; De Lucia et al. 2010). Part of the reason 
for this is the very different treatments of the winds from 
galaxies. Our extensions of the code allow us to bring the 
two approaches into closer alignment, and we include a brief 
comparison with the Galaxies-Intergalactic Medium Interac- 
tion Calculation (GIMIC Crain et al. 2009). This provides 
a series of relatively high resolution hydrodynamic simula- 
tions, featuring relatively high spatial and mass resolution 
(e = 0.5 h _1 kpc and m gas = 1.4 x 10 6 7i _1 M Q for the highest 



2 This mode is often referred to as the 'radio' mode, we prefer 
the term 'hot-halo' as it emphasises that this type of feedback is 
assumed to only be effective when the cooling time of the halo is 
sufficiently long compared to the dynamical time. See Fanidakis 
et al. (2011) for further discussion. 



resolution realisations), and that trace a representative cos- 
mological volume (four spherical volumes with comoving ra- 
dius 18/i _1 Mpc and one with comoving radius 25fe _1 Mpc). 
The simulations include radiative cooling (Wiersma, Schaye, 
& Smith 2009), and sub-grid prescriptions for star forma- 
tion and the thermodynamics of the ISM (Schaye & Dalla 
Vecchia 2008), hydrodynamically-coupled supernova-driven 
winds (Dalla Vecchia & Schaye 2008) and metal enrichment 
resulting from stellar evolution (Wiersma et al. 2009) . Some 
important successes include the X-ray scaling relations of L* 
galaxies (Crain et al. 2010), and the distribution of satellites 
and stellar halo properties of the Milky Way (Deason et al. 
2011; Font et al. 2011). Our modified GALFORM scheme 
implements similar physics, and we show that its behviour is 
very similar to GIMIC. This opens a new avenue, allowing us 
to use GALFORM to better understand how the parametri- 
sation of feedback impacts upon the formation and evolution 
of galaxies, and thus guide both the development of sub-grid 
treatments in hydrodynamical simulations, and the interpre- 
tation of observational data. 

The structure of this paper is as follows. In §2, we 
discuss the implementation of supernova-wind driven feed- 
back schemes in semi-analytic models, and introduce a new 
scheme that allows gas and metals to be expelled from low 
mass haloes and later reaccreted when the binding energy 
of the halo has increased significantly. In contrast to many 
previous models, we do not assume that material is always 
reaccreted after a certain timescale, or that expelled material 
is always lost from the hierarchy. In §3, we present a compar- 
ison of different feedback scalings, focusing on the difference 
between schemes that scale the parameters describing winds 
with halo mass, and those that adopt fixed parameters. In 
§4, we explore how feedback from supernovae can be com- 
bined with feedback from black holes in order to generate 
an exponential break in the mass function. In particular, we 
compare the effect of 'hot halo' mode feedback with that 
of strong quasar-driven winds (similar to those considered 
by Springel, Di Matteo, & Hernquist (2005)). We consider a 
number of additional observational constraints in §5. While 
several of the feedback schemes are able to reproduce the 
high mass part of the SMF, we show that the specific star 
formation rate and the downsizing of galaxy formation are 
important orthogonal constraints. We present a summary of 
our conclusions in §6. Except where otherwise noted, we as- 
sume a WMAP7 cosmology, fl = 0.272, A = 0.728, n b = 
0.0455, h = 0.704, n 3 = 0.967, cr 8 = 0.81 (Komatsu et al. 
2011). Throughout, we convert observational quantities to 
the h scaling of the theoretical model so that stellar mass 
are quoted in /i _1 M Q etc. 

2 FEEDBACK IN PHENOMENOLOGICAL 
MODELS 

2.1 Feedback as a Galactic Fountain 

We begin by reviewing the conventional approach to feed- 
back in GALFORM, focussing on the implementation used 
in Bow06. To summarise the key features, gas is expelled 
from the disk and assumed to circulate in the halo, falling 
back to the disk on roughly a dynamical time if the cooling 
time is sufficiently short. If the cooling time is long com- 
pared to the dynamical time, the halo is susceptible to a 
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hot-halo mode of feedback if a sufficiently massive central 
AGN is present. The scheme results in a galactic fountain 
with material rising from the galaxy disk and later falling 
back. In the case of haloes with short cooling times, it may 
be more appropriate to picture the circulating material as 
cool clouds rather than as material heated to the halo virial 
temperature. 

We parameterise the rate at which gas is expelled from 
the disk into the halo as 



M ou tflow = /3M» 



(1) 



where M* is the star formation rate and the macroscopic 
mass loading factor, and is 



"fdisk 
Uhot 



(2) 



Here, Vdisk is the circular speed of the galaxy disk, the pa- 
rameter i>hot sets the overall normalisation of the wind load- 
ing, and the parameter a ho t determines how the mass load- 
ing of the wind varies with the disk rotation speed. We will 
be careful to explicitly distinguish the macroscopic loading 
factor, j3, which represents the loading of the wind as it es- 
capes from the galaxy into the halo, from the sub-grid load- 
ing factor rj, used in hydrodynamical simulations to rep- 
resent the amount of interstellar medium (ISM) material 
heated or kicked by the supernova remnant. If winds are 
hydrodynamically coupled, the macroscopic mass loading is 
very likely to be significantly larger than rj. Furthermore, the 
physical processes that determine rj are themselves highly 
uncertain and their treatment in numerical models is likely 
to be resolution dependent. 

In the Bow06 implementation, material is modeled as 
leaving the disk with a specific energy comparable to the 
binding energy of the halo (ie., « W md ~ Whalo) • If we as- 
sume Wdisk ~ "halo (where Uhaio is the circular velocity of the 
halo at the virial radius), energy conservation then requires 
that ahot = 2 and that, in the standard implementation, no 
material leaves the halo completely. However, Bow06 found 
that this scaling did not sufficiently suppress the formation 
of small galaxies and a stronger scaling, ahot = 3.2, was 
adopted. This gives a good match to the observed 7f-band 
luminosity function. The stronger scaling implies that in 
small haloes supernovae couple more efficiently to the cold 
gas, resulting in a higher mass-loading of the wind. Assum- 
ing a velocity ~ i>haio/\/2 is sufficient to drive the fountain, 
the fraction of the total supernova energy needed to power 
the fountain is 



/ft 



1 / Vhat \ 

8 UsOkms" 1 / 

/ «halo 

V200 kms" 1 



(3) 



2.5 x lO^ftergMQ 1 



where esN is the energy produced by supernovae per unit 
mass of stars formed. Assuming esn = 2.5 x 10 49 h erg Mq 1 
(appropriate for a Chabrier IMF, Dalla Vecchia & Schaye 
2008) sufficient energy is, in principle, available to power 
the fountain in haloes more massive than «haio ~ 65kms~ 1 . 
Interestingly, Font et al. (2011) find that the properties of 
Milky Way satellite galaxies are best reproduced if the mass 
dependence of feedback saturates in such low mass haloes. 



2.2 Allowing for Mass Loss from the Halo 

In the revised implementation presented in this paper, we 
consider the case where the energy of the gas escaping the 
disk has systematically greater than the binding energy of 
the halo. This treatment is necessary if we are to consistently 
account for winds with outflow speeds that are independent 
of halo mass. From an observational perspective, this type 
of wind may be required in order to account for the wide 
spread distribution of metals in the universe. This type of 
wind was previously is considered in Benson et al. (2003), 
and we briefly review the implementation here. 

We introduce the parameter A cxpo i to reflect the excess 
energy of the wind relative to the binding energy of the halo. 
Specifically, we set the mean specific energy of the wind to 



E mcan — ^ AcjcDcl^halo- 



(4) 



Note that A cxpo i may be a function of halo mass (see sec- 
tion 2.3). We parameterise the fraction of material that 
is able to escape using the cumulative energy distribution, 
/e(x), where x = E eBC / E mean and E eBC is a measure of the 
wind specific energy needed to escape the halo. We will 
choose a monotonic function for Je so that /b(0) = 1, 
/e(1) = 1/e ~ 0.36 and fs — ► for large x. We will set 
the energy needed to escape the halo to E CBC — «h a i (ie., 
we assume that the escape velocity is v2«haio)- Clearly this 
is an over simplification, since the true escape velocity de- 
pends on the details of the potential, the launch radius of 
the wind, the terminal radius and the ram-pressure that the 
gas encounters. While we adopt this scaling to give a sim- 
ple interpretation the wind velocities, there is likely to be 
a systematic offset when comparing with hydrodynamical 
simulations. Since it is the ratio of -Emean and _E esc that de- 
termines the result of a model, we could rescale the wind 
speeds quoted in this paper according to the new pre-factor. 
For convenience, we can represent A cxpo i as a wind speed, 
«wind, where 



lAvind 



\ 1/2 
Aexpel ^halo 



(5) 



and we will refer to models by their wind speed; however, it 
should be remembered that this is more accurately defined 
as the specific energy of the wind, and we do not intend to 
imply that the wind necessarily has a bulk outflow velocity 
of Uwind: what really matters is the fraction of the mass of 
the outflow that escapes from the halo. Combining equations 
4 and 5, a significant fraction of the outflow will escape the 
halo if 



^wind 



> y/2 



«halo- 



(6) 



We will consider the halo mass dependence of i> w ind be- 
low, but it will be useful to normalise different models at 
a particular halo mass. For example, a fiducial halo with 
«haio = 200 km s -1 , for which A cxpc i = A expe i,200 = 1-4 cor- 
responds to wind speeds, Wwind = (200, 400, 600, 800) kms -1 . 

Material that is not expelled is added to the halo fol- 
lowing the scheme described in Bow06. This includes a de- 
lay proportional to the dynamical time before the material 
is allowed to cool again (see Bow06 for details). Material 
that escapes the halo may be later recaptured as the halo 
grows in mass. We implement this by scanning through de- 
scendant haloes in the dark matter merger tree and adding 
fE(El sc / E mcs , n ) - fsiElZ; 1 / E mc!ln ) to the reheated gas mass 
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at each step (where El sc refers to the escape energy of the 
descendant halo at timcstep i, and i ranges from the step at 
which the energy is injected to the final output time, i = 0). 
Mass that is added to the halo becomes able to cool on the 
dynamical timescale (which we define as GMhaio/Vhaio)- It 
may not be able to cool if the cooling time is long and AGN 
feedback is sufficiently effective. The step in E eBC may be 
small if the halo grows only a little by accretion, or may be 
large if the halo is accreted to become part of a much larger 
structure. 

Note that this scheme differs significantly from the su- 
perwind implementation of Baugh et al. (2005), in which 
expelled material is not considered for recapture. Since the 
overall baryon fractions of clusters of galaxies are inferred to 
be close to the cosmic abundance, recapture must be an im- 
portant part of the feedback process. Finally, we note that 
some semi-analytic models adopt a feedback scheme in which 
expelled material becomes available for cooling or star for- 
mation on a timescale that is much longer than the dynami- 
cal time. This is an approximation to the superwind scheme 
that we have described here, but it is not accurate since it 
does not take the growth rate of the halo into account. 

In order to fix on a scheme we must choose an appropri- 
ate form for the (cumulative) distribution functions /g. Ben- 
son et al. (2003) chose an exponential form, /_b = e~ x . This 
leads to a broad spread of wind particle energies. On the 
basis of their observational data, Steidel et al. (2011) sug- 
gest that wind outflow is more sharply peaked, and we also 
find that a more sharply peaked distribution better match 
the results of hydrodynamical simulations. In the following 
models we will assume Jb = exp(— x 6 ) in what follows. The 
precise choice of power is not important, however, and we 
obtain similar results for f E = exp(— x 2 ). 

2.3 A Generalized Feedback Model 

In contrast with GALFORM, and largely because of a lack 
of observational motivation for any particular scaling, hy- 
drodynamical simulations have mostly adopted the simplest 
case of assuming that the mass loading and velocity of winds 
are independent of halo properties. We can easily adapt 
the revised feedback implementation described above to in- 
vestigate such a scheme in GALFORM, by explicitly in- 
cluding a halo-mass dependence in Eq. 4. We will adopt 
fhaio = 200 kms -1 as a fiducial halo mass at which to com- 
pare the wind mass loading and wind speed for models with 
different a cxp so that we express the mass dependence of 
the wind speed as A cxpc i = A exp ei,200 (u lm io/200kms~ 1 ) Q:oxp , 
where a oxp is a dimensionless parameter that differentiates 
different feedback models. Combining this with Eq. 4 and 
Eq. 5 gives 

Qioxp+2 

(7) 



22/ 

^wind = «wind,200 ^7 



Uhalo 



,200 kms- 1 
For the mass loading we have 

a a ( «wind,200 \ ""hot 

^^'Uokms-J (8) 
In the case a oxp = we recover the wind specific energy 
scaling with the specific binding energy of the halo. In the 
case a cxp = —2, the wind speed is independent of halo mass. 
Unless otherwise stated, we will set 



so that a fixed fraction of the supernova energy is used to 
drive winds in haloes of all masses (assuming Udisk ~ «haio 
in Eq. 2). 

Given a wind speed and mass loading normalisation, 
v wind, 200 and /?200, the code parameters are 



A 



e X pel,200 



lVind,200 

200 kms" 1 



and 



fhot 



/?20o 1/ahot 200kms" 



(10) 



(11) 



This allows simple comparison to older models. Note that 
the original GALFORM parameter, Vhot, expresses the mass 
loading, and is not a measure of the specific energy of the 
wind. To avoid this confusion in this paper, we will use the 
macroscopic mass loading parameter P200 to label models in 
what follows. The maximum available supernova energy sets 
a limit /3 2 oo(Vwind,2oo/200 kms -1 ) 2 < 32. 

If a cxp = —2 (ahot = 0), the macroscopic mass loading, 
/3, is independent of the halo potential. This would be the 
case if the wind were completely decoupled until it that had 
escaped from the halo. This would make it impossible to 
frame the feedback in terms of the standard GALFORM 
parameters, and for this paper, we will assume that there 
is always a small coupling between the wind loading and 
the halo mass. We adopt a cxp = —1.9 (a no t = 0.1) as our 
minimum value. 



While a c 



and —2 are natural choices, there is 



Ohot 



+ 2 



(9) 



no a- priori reason to adopt a particular value of a exp , and 
we will consider a cxp = — 1 («hot = 1) as an intermediate 
value. For these parameters, the speed of the wind scales 
with y/Vhsdo, and its mass loading scales as 1/vhaio- In this 
case, the material expelled from smaller galaxies is more 
likely to escape the halo, but this is a weaker function of 
mass than in the superwind case discussed above. The scal- 
ing of the wind mass loading is similar to the momentum 
driven model used by Oppenheimer & Dave (2008), but note 
that we will assume that the ratio of the total energy of the 
wind (not its total momentum) to the mass of stars formed 
is independent of halo mass. 



2.4 Parameter values 

The parameters of the best fitting models are given in Ta- 
ble 1. To simplify comparison with previous work, we have 
translated the feedback parameterisation in Bow06 (fhot = 
485 kms -1 , a ho t = 3.2, A expe i,200 = and a coo i = 0.58) 
into to the more generalized parameters considered in this 
paper. Where parameter values are not explicitly given, we 
adopt those in Bow06 with the exceptions given below. First, 
we now adopt a background cosmology that is consistent 
with the WMAP 7-year results (Komatsu et al. 2011). Sec- 
ondly, we use a stellar yield of p y i c id = 0.04 in order to 
improve the match of galaxy colours as discussed in Font et 
al. (2008), and a default halo gas distribution with core ra- 
dius of 0.025r v i r as discussed in Bow08. With these revisions 
we make small shifts in the standard feedback parameters 
in order to restore a good fit to the local if-band luminosity 
function. 

The parameters of the base-line model are given in 
the second row of Table 1. We use v bot = 430 kms -1 and 
Qcooi = 0.52, where Q CO oi determines the ratio of free-fall and 
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Model 




"hot 


/3200 


?Avmd,20o/kms 1 


f^cxp 


^cool 


Bow06 


Bow06 


3.2 


17 






0.58 


+ W7 cosmology 


B8W7 


3.2 


12 




0.0 


0.52 


Superwind 


pGIMIC 


0.1 


8 


275 


-1.9 






SW 


0.1 


8 


180 


-1.9 


(0.35) 


Momentum Scaling 


MS 


1.0 


8 


200 


-1.0 


(0.45) 


Energy Scaling 


ES 


2.0 






0.0 





Table 1. The table shows the models considered in this paper, and compares to the values adopted in Bow06. Up dating the comsological 
paramaters requires that we make small adjustments to restore a good match to the observed mass function. We will refer to this model 
as B8W7. The following models correspond to the optimal paramaters for the different feedback schemes considered in the text. The 
values in brackets indicate the values adopted for the AGN "hot halo" of feedback in §4.1. 



cooling times at which haloes are taken to be hydrostatic (as 
opposed to being classified as 'rapidly cooling') so that only 
when t C ooi(Vcooi) > a~ 1 ol M r c°°i) is the AGN feedback effec- 
tive. These differences have little impact on the properties 
of sub-L* galaxies. Since we are initially concerned with the 
faint end of the luminosity function, we begin by disabling 
the AGN feedback scheme by setting a coo i = 0. This allows 
us to make a simple comparison to hydrodynamical calcula- 
tions that do not include AGN feedback. 

We consider three supernova-driven feedback models 
following the feedback scheme discussed in the previous sec- 
tion. In what follows, we refer to these as 'superwind' (SW), 
'momentum scaling' (MS) and 'energy scaling' (ES) mod- 
els. This nomenclature reflects how the mass loading scales 
with halo mass, corresponding to ahot = 0.1, 1 and 2. For 
each case, we typically consider six values of the mass load- 
ing at i>hak> = 200 kms -1 , /32oo = 1,2,4,8,16 and 32, and 
five values of the wind launch velocity at this halo mass, 
«wind,20o = 100, 200, 400, 600 and 800 kms -1 . For the super- 
wind and momentum scaling models, optimal values of the 
feedback parameters have been chosen to provide a reason- 
able match to the SMF above M» ~ 10 9 ft, _1 M Q (a par- 
ticularly good fit was not possible for the ES scaling) and 
the corresponding parameters are given in Table 1. We also 
consider a model intended to replicate the supernova-driven 
winds implemented in the GIMIC hydrodynamical simu- 
lations (Crain et al. 2009). These adopt a fixed sub-grid 
wind mass loading of r\ = 4, and a launch wind velocity of 
600 km s" 1 . This requires ~ 80 percent of the available su- 
pernova energy being used to drive the wind. We find that 
we can best reproduce the resulting stellar mass function 
produced by the simulation if we adopt a somewhat lower 
macroscopic wind speed, u W md,200 = 275kms _1 , and higher 
mass loading, P200 = 8 (although a broad range of parame- 
ters give similar results). Since the winds in the simulation 
are always hydrodynamically coupled, the increased mass 
loading of the wind is not surprising. We will refer to this 
'pseudo-GIMIC' model as pGIMIC below. 

While we initially consider the models in the absence 
of AGN feedback, we will later adjust the AGN feedback 
parameter a coo i in order to set the mass scale of the break 
in the luminosity function. For this we use the AGN feedback 
scheme of Bower, McCarthy, & Benson (2008). This allows 
material to be expelled from the central regions of haloes 
through AGN heating, and results in the baryon fraction 
in galaxy groups being much lower than the cosmic average. 
This provides a much improved match to the observed X-ray 
properties of these systems. The values of a coo i that give a 




9.5 10,0 10.5 11.0 11.5 
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Figure 1. Comparison of the stellar mass function of the Bow06 
model (blue line) with the base-line B8W7 model used in this 
paper (green line). This is based on the WMAP7 cosmology and 
includes AGN "hot halo" feedback following Bower et al. 2008. 
The two models are almost indistinguishable. To illustrate the 
importance of AGN feedback, we show the effect of turning off the 
AGN feedback (red line). We also show the effect of then adopting 
the «hot = — 2 (cyan line). For comparison, observational data are 
shown as black points with error bars. The data is taken from Bell 
et al. (2003) (circles) and Li & White (2009) (crosses). 



good match to the observed luminosity function are given in 
brackets in the table. The ES model could not be adjusted 
to give a sufficiently flat faint end slope to the SMF, and we 
do not consider the role of AGN in this model. 



3 GALACTIC WIND MODELS 

3.1 Conventional GALFORM Winds 

We use the galaxy SMF as the starting point for compar- 
ing the models and the data. We use the determination of 
the stellar mass function by Bell et al. (2003) and Li & 
White (2009), correcting the IMF to the Kennicut param- 
eterisation. We convert the observational data to the h~ x 
dependency on the Hubble parameter of the the theoreti- 
cal models. The data are compared with the Bow06 model, 
and the baseline B8W7 model in Fig. 1. Both models as- 
sume that supernova-driven winds do not escape the parent 
halo; instead the material ejected from the disk circulates 
in the halo and returns on the dynamical timescale if the 
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cooling time is short. The two models are almost indistin- 
guishable, since changes in the background cosmology and 
the AGN feedback implementation have been intentionally 
compensated for by small changes in the feedback param- 
eters. As expected, both provide a good description of the 
observational data. 

When we consider alternative wind descriptions, we will 
not initially consider AGN feedback. To establish a baseline 
for the comparison in the absence of AGN feedback, we also 
show the B8W7 model with AGN feedback disabled (by set- 
ting a C ooi = 0). This is shown as a red line in the figure. 
The roll-over of the galaxy luminosity function is almost 
non-existent in this model, while the SMF is largely unaf- 
fected below 10 11 M Q . This is encouraging, since it shows 
that the two processes involved in matching the shape of 
the galaxy luminosity function (ie., eliminating the over- 
abundance of galaxies of fainter than ~ 10 h~ Mq and 
reducing the abundance of galaxies above the break in the 
mass function) can be separated. We therefore focus our ini- 
tial discussion on the modes of supernova-driven feedback, 
and consider models that do not include an AGN feedback 
component. 

The Bow06 and B8W7 models achieve a good fit to the 
abundance of low mass galaxies because of the very strong 
mass dependence of the feedback (a hot — 3.2). For compari- 
son, we show a model with feedback parameters that might 
be considered to be better motivated by theory. By selecting 
«hot = 2, the wind speed is tuned to the binding energy of 
the halo. As can be seen, this results in a somewhat steeper 
faint end-slope and a relatively poor match to the observed 
stellar mass function. In the following section we consider 
in more detail how the choice of feedback scheme affects the 
mass function, and how the wind parameters can be opti- 
mised to improve the fit. 

Although the SMF provides a good way to compare 
the results of different feedback schemes with observation, 
it is far from simple to interpret the changes in terms of 
the effect of the different feedback schemes. For example, 
increasing the effectiveness of the feedback scheme shifts 
galaxies to lower stellar mass, and so only affects the nor- 
malisation of the mass function indirectly. The suppression 
of the normalisation arises both because of the lower abun- 
dance of the haloes of greater mass and because of the range 
of halo masses that contribute galaxies of a particular stellar 
mass. A better way to understand the effect of the schemes 
is therefore to plot the stellar mass of the central galaxy as 
a function of halo mass. Since the scatter is not strongly 
constrained observationally (Moster et al. 2010), we use the 
relationship found in the Bow06 and B8W7 models as a best 
guide to the relationship expected in the real Universe. We 
can then understand how various feedback schemes affect 
this relationship, and compare with constraints yielded by 
abundance matching observations with theoretical subhalo 
mass functions. 3 Of course, the parameters chosen in Bow06 
and B8W7 are not unique, and other parameter combina- 
tions can give similar quality fits to the mass function and 
other datasets, as we have shown in Bower et al. (2010), but 
the models to provide a well documented starting point for 
our comparison of different feedback schemes. 

The first panel of Fig. 2 illustrates the dependency of 
central galaxy stellar mass on halo mass for the Bow06 
model. The red line shows the mean relation. The dashed 



black line shows a broken power law approximation to the 
model, described by 

logM, = 10.1 + 0.54(log(M ha io) - 11.7) 

if log(Mhalo) > H.7 

= 10.1 + 2(log(M ha io) - 11-7) 

iflog(M h aio) < 11-7. (12) 

At the break in the curve, 14 percent of the baryons in the 
halo have been converted into stars. Note that there is con- 
siderable scatter about these relations in the model. We re- 
peat this relation in all plots so that they can be compared 
easily. The error bars indicate the ±1<7 range of the model 
galaxies. Close to the break, the scatter in the model exceeds 
an order of magnitude. We supplement the empirical approx- 
imation with relations from Moster et al. (2010) and Guo et 
al. (2010). The relations shown are derived from matching 
the abundance of sub-haloes in TV-body simulations to ob- 
servational data, assuming that the scatter in the relation is 
negligible. The models are based on (fi m , Qa, h, as, n s ) — 
(0.26,0.74,0.72,0.77,0.95) and (0.25,0.75,0.73,0.9,1) cos- 
mologies respectively, but note that the differences in the 
predicted abundance of 10 10 — 10 12 /i _1 M© haloes are small. 
Thus the relations are similar for stellar masses below 
lO^r'Me but are offset from the Bow06 relation at high 
masses due to the large scatter about the mean relation. Be- 
cause of the steep break in the mass function, scatter boosts 
the abundance of massive galaxies relative to a relation with- 
out scatter (see Moster et al. (2010) for further discussion), 
and thus the scatter and the normalisation of the high mass 
M* — -Mhaio are tightly correlated. 

The second panel of Fig. 2 shows the baseline B8W7 
model. The scatter in M* at a given halo mass is reduced 
compared to Bow06, although it is still larger (a = 0.3 dex) 
than the scatter explored by Moster et al. (2010) (up to 
0.15 dex), particularly around the break in the relation. The 
final panel shows the effect of disabling AGN feedback in the 
B8W7 model. The power-law relation now extends to higher 
mass before slowly rolling over as the result of the increasing 
cooling times of massive haloes. The scatter in the relation 
around Mhaio ~ 10 12 ft _1 Mo is now much reduced. This 
arises because the efficacy of AGN feedback in this model 
has a strong dependence on the accretion history of haloes 
(sec Bow08 for further discussion). 

In summary, the B8W7 model provides a match to the 
observed stellar mass function due to the very strong halo 
mass dependence of the wind mass loading and the suppres- 
sion of cooling in haloes with relatively long cooling times. 
This paper investigates whether models with more general 
feedback schemes can achieve a similar success. 



3.2 Superwind Models with Fixed Wind Speed 
and Mass Loading 

3. 2. 1 Effect of wind parameters 

We now consider superwind (SW) models, in which the mass 
loading and velocity of winds are (almost) independent of 
the halo mass. This mimics the schemes that are commonly 
adopted in hydrodynamical simulations. We begin by con- 
trasting the results with the baseline B8W7 model. 

Fig. 3 shows the effect of varying the mass loading for 
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Figure 2. The M„ - Mh a i relation for central galaxies. In this figure the panels show (from left to right) Bow06; the base-line model 
B8W7; and B8W7 without AGN feedback. Blue points show a random selection of haloes from the model, while the solid red line shows 
the median, and error bars show the la scatter. The solid black line shows an empirical fit to the Bow06 model, while black dashed 
and dot-dashed lines show halo abundance matching models from Moster et al. (2010) and Guo et al. (2010) respectively. The models 
in the first two panels include AGN feedback and successfully match the observed stellar mass function. These are well characterised 
by a broken power-law. Although the first two panels have slightly different background cosmology, the relation between halo mass and 
stellar mass is similar. Disabling AGN feedback, as shown in the last panel, results in a much weaker roll-over of the relationship. 




0.0 0.5 10.0 10.5 11.0 

log^tellarMass/ft- 1 ^;,) 



Figure 3. Effect of mass loading in the superwind (SW) model. 
These models (solid lines) have a wind velocity of w w ind 200 = 
400kms — 1 , and mass loading /?200 = 1,2,4,8,16. We also in- 
clude a model with low mass loading (/?200 = 1) and low speed 
(''wind 200 = 100 kms -1 ) to show the effect of suppressing the 
feedback channel (blue line). With low mass loading, the faint 
end of the function is extremely steep. Increasing the mass load- 
ing suppresses the normalisation of the low mass part of the mass 
function, but tends to introduce a noticeable dip in the SMF if 
the loading is high. With a high mass loading also suppresses the 
abundance of the highest mass galaxies, but shape of the mass 
function does not match the observational data. Since AGN feed- 
back has not been included, the cut-off seen here is the result of 
the long cooling times in larger haloes. 



a fiducial wind speed of u w ind,200 = 400 km s -1 . For com- 
parison, the blue curve shows the effect of low mass loading 
and low wind speed such that supernova driven feedback 
is ineffective. Increasing the mass loading reduces the nor- 
malisation of the mass function below M*, but the power- 



law dependence at low galaxy luminosities becomes steep, 
increasing the discrepancy with observations. As the mass 
loading increases above P200 = 4, feedback takes a 'bite' 
out of the mass function. This can be understood as a tran- 
sition in the effectiveness of feedback. In high mass haloes 
(fhaio > «wind,20o) material falls back onto the central galaxy 
on the dynamical timescale, while in low mass haloes it is 
expelled and ceases to be available to fuel star formation. 
The timescale for the return of expelled material therefore 
makes a transition when the two speeds are equal (see Op- 
penheimer & Dave 2008). (The effect is most clearly seen by 
plotting the stellar mass of galaxies against their halo mass, 
as we discuss below.) Although increasing the mass loading 
tends to suppress the abundance of M* galaxies, the faint 
end slope of the mass function is always much steeper than 
the observations. It is not possible to improve the fit to the 
mass function by adjusting this parameter. 

Although the abundance of bright galaxies greatly ex- 
ceeds the observations, a roll-over is evident when the high 
mass loading. As was the case where AGN were disabled in 
the B8W7 model, the roll over is driven by the long cooling 
times of high mass haloes. In this situation, the bottle-neck 
is the cooling time of the material in the halo, and the star 
formation rate is proportional to the inverse of the wind 
mass loading. 

The effect of changing the wind speed is explored in 
Fig. 4, where we show the mass function obtained by vary- 
ing v w ind,200 = 100 — 400kms _1 at a fixed mass loading 
of /?200 = 8. At a wind speeds greater than 200 km s -1 , 
the mass function begins to resemble the observational data 
with a fiat slope around the knee of the luminosity function. 
This is encouraging: we infer that that a suitable choice of 
wind parameters enables a match to the properties of 10 9 ' 5 
- 10 11 /i -1 Mq galaxies to be obtained. However, while the 
model matches the observational data down to M, ~ 10 9 ' 5 , 
the number of galaxies rises rapidly at lower masses, and 
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Figure 4. The effect of changing the wind speed in the superwind 
(SW) model. These models have a fixed mass loading /3200 = 8 
and varying wind speed u w i n d = 100 — 400 km s -1 . The figure 
shows how changing the speed of the wind adjusts the location 
and depth of the dip in the SMF, although the faintest and most 
massive galaxies are almost unaffected. An intermediate wind 
speed has the effect of producing a flat portion of the luminosity 
function. 

an additional feedback mechanism would need to be intro- 
duced to explain the low abundance of dwarf galaxies. Since 
the halo masses of these galaxies are sufficiently high that 
they are unlikely to be affected by photo-heating (Crain et 
al. 2007; Okamoto, Gao, & Theuns 2008), the only option 
would be to explore winds that scale with halo mass. On 
the other hand a high abundance of faint galaxies (at z > 6) 
would provide an abundant source of ionising photons to 
drive the re-ionisation of the universe (Benson et al. 2006; 
Jaacks et al. 2011). If the mass-dependent scheme suppress 
the formation of pre-reionisation small galaxies too dramat- 
ically there will not be sufficient photons to re-ionise the 
universe. The constraint is quite weak, however. Even with 
the strong halo mass dependence winds in the B06 model, we 
find that it is sufficient to assume that feedback saturates 
when «h a i < 65kms _1 in order to provide the necessary 
ionising flux (Font et al. 2011). 

In summary, with suitable choice of parameters, the SW 
scheme offers an attractive explanation for the flat portion of 
the SMF in the range 10 9 5 - 10 11 /i _1 Mq. The model, how- 
ever, predicts that lower mass galaxies will be more abun- 
dant than observed. In contrast, the strong B8W7 halo mass 
dependence of feedback in the B8W7 model results in a flat 
stellar mass function down to below 10 s fc _1 Mg. 

3.2.2 Comparison with Hydrodynamical Simulations 

The SW feedback scheme is similar to the approach adopted 
in many hydrodynamic simulations, and it is interesting 
to briefly compare the results in order to gain confidence 
that the phenomenological description we use appropri- 
ately represents the physics of a full hydrodynamical treat- 
ment. We first compare with the Galaxies and the Inter- 
galactic Medium Calculation (GIMIC, Crain et al. (2009)). 
The highest resolution realisations of these simulations have 
gas particles with mass 1.45 x 10 6 /i -1 Mq, and a softening 



Figure 5. Comparison with the SMF of the GIMIC hydro- 
dynamic simulation (dashed line) and the GALFORM model 
(pGIMIC) with (« w i n d,200, A200) = (275,8) (solid blue his- 
togram). The GALFORM parameters have been chosen to re- 
produce the position and amplitude of the dip seen in the hydro- 
dynamic simulation. The GIMIC wind parameters are (600,4), 
but hydrodynamical coupling is expected to result in additional 
macroscopic mass loading and lower effective wind speeds. Ob- 
servational data is reproduced from Fig. f . 

length of 0.5/i _1 kpc. This is sufficient to resolve the on- 
set of the Jeans' instability in galactic disks while at the 
same time allowing reconstruction of a representative cos- 
mological volume. Feedback is implemented by imparting 
kinetic energy to stochastically chosen neighboring particles 
of newly-formed stars. The kicked particles remain hydrody- 
namically coupled at all times. The simulations adopt sub- 
grid wind parameters r\ = 4 and « w ind,200 = 600kms _1 . 
This choice was motivated by observations suggesting that 
the wind speed was independent of halo mass (see Martin 
(2005) for discussion) and by requiring that the total stel- 
lar mass density matched the observed universe. Alterna- 
tive choices are explored in the OWLS simulations (Schaye 
et al. 2010). These parameters determine the input proper- 
ties of particles. Since the simulation is fully hydrodynamic 
(so that wind particles remain hydrodynamically coupled to 
the surrounding gas particles), we should not expected them 
to directly translate into the macroscopic wind parameters 
used in GALFORM. 

In Fig. 5 we compare the SMF from GIMIC with the 
SW model (a full comparison of individual galaxy merger 
trees will be presented in a future paper). It is important 
to note that the GIMIC simulations did not include AGN 
feedback, and were constrained to matching the observed 
stellar mass density, rather than the portion of the mass 
function below stellar masses of 10 11 h _1 M . In order to 
run the GALFORM model, we revert to the cosmologi- 
cal parameters used in Bow06 (both models are based on 
the Millennium simulations, Springel, Di Matteo, & Hern- 
quist (2005)). With suitable choice of wind normalisation 
(iWd,2oo = 275 kms" 1 and /3 2 oo = 8), the GALFORM code 
reproduces the hydrodynamic mass function well. Although 
the GALFORM wind has lower speed and somewhat higher 
mass loading, it should be remembered that these are the 
effective macroscopic wind parameters. As was shown by 
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Dalla Vecchia & Schaye (2008), the ram pressure induced 
by the hydrodynamic coupling of winds tends to slow the 
outflow and increase its mass loading as it leaves the disk of 
the galaxy. 

Oppenheimer & Dave (2008) (hereafter Op08) also 
present similar models with which we can compare. These 
simulations have significantly lower mass resolution than 
GIMIC, with gas particle masses of up to 1.5 x 10 7 /i _1 Mq 
and a softening length of 1.9 ft _1 kpc. Consequently, these 
hydrodynamic simulations do not attempt to resolve physics 
within galaxies. Star formation is implemented following the 
sub-grid scheme of Springel & Hernquist (2003): winds are 
implemented kinetically, but the affected particles are decou- 
pled from hydrodynamic forces until the surrounding density 
is low. 

Op08 consider three distinct models. A high wind speed 
model, a model with much lower wind speed, and another in 
which the wind speed (and mass loading) scale (inversely) 
with the local halo velocity dispersion. Their strong wind 
model (f w ind,20o, /3'200 = 680,2) produces results that are 
similar to those of the GIMIC simulations. However, the 
'slow wind' model (w w ind,200, P200 = 340,2) provides a good 
match to the observed mass function over the range plot- 
ted in their paper. There are three regimes of the SMF 
for the slow wind. A steep slope at low mass, flat around 
10 10 fe _1 M0 and then steep again at higher masses. As we 
have seen, GALFORM can reproduce this behaviour if the 
wind velocity is low (200 kms -1 ) and the mass loading some- 
what higher (between 4 and 8). These have comparable to- 
tal wind energy to the Op08 models. In low mass haloes, 
even the high wind loading considered does not sufficiently 
suppress star formation, compared to the observations. At 
intermediate masses, the slope is roughly flat as the wind 
becomes less effective and eventually stalls. Then at very 
high mass, cooling becomes inefficient and the slope steep- 
ens. Obviously, as with GIMIC, the match to the observed 
SMF at such high masses is poor because the simulations 
do not include AGN feedback (but see Gabor et al. (2011)). 
We will consider models in which the wind parameters scale 
with the properties of the halo in Section 3.3. 

In summary, this brief comparison shows that the ex- 
pulsion scheme implemented in our phenomenological model 
describes the effects seen in hydrodynamic simulations well. 
By using these models to better explore the parameter 
space of galaxy feedback, we can create a closer connec- 
tion between phenomenological models and fully hydro- 
dynamic simulations. Comparison with the observational 
data highlights two important issues: firstly, the steep slope 
of the mass function below M„ = 10 10 ft _1 Mo and sec- 
ondly the over-abundance of galaxies more massive than 
M* = 10 11 ' 5 h~ 1 M Q . In the following sections, we will ex- 
plore how these discrepancies can be resolved by introducing 
more feedback schemes that scale with halo mass, and in- 
cluding feedback from AGN. 

3.2.3 Stellar Mass as a Function of Halo Mass 

In order to better understand how the feedback scheme can 
shape the SMF, it is useful to examine the relation between 
the halo mass and the stellar mass of central galaxies. We 
presented the relation for the Bow06 model in Fig. 2 and 
showed that the relationship can be characterised by a bro- 



ken power law. Fig. 6 illustrates the effect of changing the 
feedback scheme to the SW model. 

The first panel shows the effect of including only min- 
imal feedback, (f W md,2oo, ftoo) = (100,1). The discrepan- 
cies compared to the empirical relationship (black lines) are 
evident. At all halo masses, the associated stellar mass is 
too high, and the relation shows little change of slope. This 
model corresponds to the weak feedback (blue) line in Fig. 3. 
At a given stellar mass, the galaxies are formed in lower mass 
haloes than indicated by the observed relation. These haloes 
are much more abundant, and thus the SMF is normalised 
too high. The slope of the relation is also in clear disagree- 
ment and this results in the overly steep faint end slope of 
the predicted mass function. 

The remaining panels illustrate the effect of increasing 
the wind speed at a fixed mass loading. In the second panel, 
(«wind,2oo, P200) = (100,8). This model is shown as the blue 
line in Fig 4. The outflow has a low speed, so little mass es- 
capes the halo, but the high mass loading results in effective 
suppression of galaxy stellar mass. As a result, the model 
matches the normalisation of the knee of the SMF well, and 
this is reflected by the M„ — Mhaio relation coming close to 
the kink of the observed relationship. However, several dis- 
crepancies from the observed relationship remain clear. In 
particular, the relation shows little change of slope. While 
the relationship at high stellar mass can be improved with 
AGN feedback, the relation at lower stellar masses is too 
shallow. As a result, galaxies of a given stellar mass are over 
abundant compared to the observed relation, as is evident 
in Fig 4. 

The third and fourth panels show the effect of increas- 
ing the mass loading in the model, (-u W md,200, P200) = (200, 8) 
and (400, 8). These are shown as the green and purple lines 
in Fig. 4. The increasing the wind speed creates a kink 
in the M*-Mh a io relation, with the stellar mass formed in 
haloes around 10 11 /i _1 M Q being very strongly suppressed. 
In the kinked region, the steepening of the M*-Mhaio re- 
lation means that a particular halo mass contributes to a 
wide spread of stellar masses, resulting in a suppression of 
the mass function normalisation. By suitable adjustment of 
the parameters, the suppression can be tuned to create a flat 
portion of the SMF. We will exploit this in Section 4.1. Be- 
low the kink, however, the slope of the M*-Mh a io relation is 
much shallower than that seen in the B8W7 model (Fig. 1) 
and the slope of the mass function is therefore inevitably 
steeper than the observational data. 

The kink is created by the wind stalling at a particu- 
lar halo mass so that material no longer escapes the halo 
and falls back on the dynamical timescale. Increasing the 
speed of the wind shifts the region of the kink, but leaves 
the relations at high and low halo masses unchanged. This is 
reflected in the SMF, with the wind speed effecting a transi- 
tion region between unchanging abundances of high and low 
mass galaxies. The transition between the regimes appears 
to become steeper for higher wind energies, resulting in a 
noticeable dip in the mass function. 

3.3 Momentum Scaling Models 

While Op08 find that the slow wind model fits the observed 
data over a range from 10 9 ' 5 to 10 11 7i -1 M , their preferred 
model is one in which the feedback parameters vary with 
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Figure 6. Comparison of the M* - Afh a i relation for SW models with various parameters. The first panel shows a model with weak 
feedback (w w i n d,200i P200) = (100, 1). The remaining panels show (i> w ind,200i 02oo) = (100,8), (200,8) and (400,8) (from left to right) to 
illustrate the effect of increasing the wind speed at a fixed mass loading. The SMFs corresponding to these models can be seen in Fig. 4, 
the effects seen in the SMF are more readily interpreted in terms of the amplitude and slope of the M, — M^aio relation. Black lines 
show the expected relationship derived from on observational data (see Fig. 2). 
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Figure 7. The effect of using the "momentum" scaling (MS) model (cthot = 1)- The first panel shows a wind speed normalisation of 
"wind = 200 km s _1 , with varying mass loading, /3200 = 2 — 32. For a low wind speed, a high normalisation of the mass loading effectively 
suppresses the formation of low mass galaxies leading to an encouraging match to the a large portion of mass function. The second panel 
shows a fixed wind loading of 8 and varying normalisation of the wind speed, ^ w ind = 100 — 400kms — 1 . Increasing the wind speed has 
the undesirable effect of creating a dip in the mass function. 



halo mass. Their preferred scheme is intended to mimic the 
effect of a momentum driven wind (Murray, Quataert, & 
Thompson 2005), and the ratio of total momentum to mass 
of stars formed is held fixed. We adopt a similar scheme, scal- 
ing the wind mass loading inversely with the circular speed 
of the disk (ie., /3 oc 1/tidisk)- In contrast to Op08, however, 
we scale the wind speed so that the total wind energy (per 
stellar mass formed) is independent of halo mass. (The to- 
tal of the wind momentum in Op08 exceeds the momentum 
available from photon by almost an order of magnitude). 
The effect of using this 'momentum scaling' (ai lot = 1) in 
the GALFORM model is shown in Fig. 7. The first panel 
shows a model with relatively modest wind speed normali- 
sation (u w ind,200 = 200 kms" 1 ), considering a range of mass 
loading normalisations — 2 — 32. Note that in low 



mass haloes, the mass loading will be higher and the wind 
speed lower. The panel shows that normalisation of the SMF 
steadily decreases as the wind speed increases. The overall 
SMF is flatter than that seen in the SW case (where the 
wind properties are independent of halo mass, Fig. 4), and 
we see that some models compare favourably with the data 
above a stellar mass of 10 9 /i^'Mq. 

The second panel shows the effect of increasing the wind 
speed at a fixed mass loading, /?200 = 8. As the wind speed 
increases, the abundance of Af» galaxies is suppressed. As 
we have seen in the SW models, a high mass loading can 
create a dip in the mass function (where u w ind ~ Vhalo)- The 
dip tends to be more smeared out for the MS wind, how- 
ever. The origin of this feature is shown clearly in Fig 8, 
where the M» - Mhaio relation is shown for an MS mod- 
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Figure 8. Comparison of the M* - Afhalo relation for two 
momentum scaling models (oihot = !)• The first panel shows 
(^wind^oOi A200) = (150,8), the second panel (400,8). Compar- 
ing this plot with Fig. 6 shows that the principle effect of the 
halo mass dependence of the wind is to smooth out the kink in 
the M„ - Afhalo relation. See Fig. 2 for explanation of lines and 
symbols. 



els with parameters (u w ind,200, P200) = (150,8) and (400,8). 
The effect of the MS feedback scheme is to introduce a kink 
into this relation, with the location of the kink depending 
on the wind speed normalisation. The effect is similar to 
that seen previously in the SW models, but the kink is more 
diffuse, resulting in a smoother transition between the high 
and low mass regimes. The second panel of this figure should 
be compared with the last panel of Fig 6 as the feedback is 
the same in 200 km s _1 haloes in both cases. The kink in the 
M* - Mhaio relation is clearly much smoother in momentum 
driven model, resulting in a less prominent dip in the mass 
function. Below M» = 10 9 h~ 1 M Q , the slope of the M« - 
Mhaio relation is again relatively shallow leading to an over 
abundance of low mass galaxies. 

The difference between the MS and SW models is also 
seen in the behaviour of the high mass end of the SMF. At 
high masses, little of the feedback material is able to escape 
the halo, but the mass loading of the two models differ. As 
a result, for models with equal normalisation, the stellar 
mass associated with high mass haloes in higher in the MS 
model than in the SW model (compare the last panels of 
Fig. 8 and 6, for example). This is reflected in a greater 
abundance of high mass galaxies in the MS model vs. SW 
(compare Fig. 7 and 4). We will show, later, that this has 
important consequences for the abundance of the massive 
galaxies at high redshift. 

In summary, the dependencies on mass loading and 
wind speed in the MS model show similar trends to the SW 
model. However, the rise in the abundance of the faintest 
galaxies is shallower, and the dip in the SMF tends to be 
smoothed out. With suitable choice of wind parameters, this 
model is able to match the observed stellar mass function 
over a greater range of galaxy mass. 



Figure 10. Comparison of M, vs. M^aio for a model with energy 
scaling of the wind (ES, «hot = 2). The two panels show models 
with (i>wind,200>ft!Oo) = (200,1) and (200,32). These span the 
range of models shown in the right hand panel of Fig. 9. See 
Fig. 2 for explanation of lines and symbols. 



3.4 Energy Scaling Models 

Finally, we consider models in which the wind speed is a 
fixed multiple of the halo circular velocity, such that a fixed 
fraction of the wind escapes regardless of the halo mass. In 
the left hand panel of Fig. 9, we fix P200 = 4 and allow 
the velocity of the wind to increase, w W md,20o = 100 - 400. 
The results for normalised wind speeds between 100 and 
200 (blue and green lines) are identical because little mate- 
rial has sufficient specific energy to escape the halo. Further 
increases in wind speed change the normalisation dramati- 
cally as material leaves the halo and takes longer to become 
available for cooling again. However, the shape of the SMF 
changes little, and it is not possible to recreate the dip in the 
mass function that was seen in previous models. In the right 
hand panel, we show the effect of varying the mass load- 
ing of the wind. For this model, the effect is similar to that 
of varying the wind speed. Since there is no characteristic 
mass at which the wind stalls, the loss of material from the 
halo is similar regardless of whether a relatively small frac- 
tion of baryons are expelled with high specific energy (and 
thus remain outside of the halo for an extended period), or 
a large fraction of material is expelled with lower specific 
energy. Finally, we note that the abundance of high mass 
galaxies continues the increasing trend seen Figs. 7 and 4 
due to the decrease in mass loading in high mass haloes (for 
equal normalisation). 

Figure 10 shows the effect of this type of feedback on 
the M*-Mhaio relation. As can be seen, the kink that en- 
abled us to produce a flat component of the mass function 
in the first two feedback schemes is absent. The scaling of 
the wind speed with halo mass means that winds do not stall 
at a particular halo mass. Instead, the feedback steepens the 
overall slope of the M*-Mh a io relation. However, although 
it is closer to the observed relation, the difference in slope 
leads to a significant mismatch with the SMF, as shown in 
Fig 9. This clearly illustrates the way in which the normal- 
isation of the mass function is strongly dependent on the 
slope of the M*-Mh a io relation as well as its normalisation. 
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Figure 9. (left panel): The effect of varying the wind speed in the energy scaling (ES) model («hot = 2). These models have increasing 
wind speed normalisation, v w i n d 200 = 100, 200, 300, 400; all the models shown have @200 = 4. Increasing the wind speed normalisation 
results in further suppression of the SMF. However, the effect is strongest around the knee of the mass function, and the resulting function 
does not match the observational data, (right panel): The effect of varying the wind loading at a fixed wind speed. For this model, the 
effects of varying mass loading and wind speed are quite similar and the mass function cannot be shaped to match the observational 
data. 



Overall, the effect of this type of feedback is less en- 
couraging, and we do not consider this model further. The 
effects of mass loading and wind speed are very similar, and 
the primary effect of both is to suppress the normalisation 
of the stellar mass function rather than to alter its shape. 
While including AGN feedback induces a break at the bright 
end of the mass function by suppressing cooling in hydro- 
static haloes, the faint end slope is not affected. In contrast, 
the B8W7 model achieves a much improved match to the 
mass function by adopting a stronger halo mass dependence 
of the wind mass loading. 



4 THE ROLE OF AGN FEEDBACK 

4.1 The hot-halo (or "radio") mode 

None of the models discussed so far are able to match the 
abrupt turn over of the SMF. In this section, we consider 
the "hot-halo" mode of AGN feedback, associated with the 
heating of the group and cluster diffuse material. Bow06 ar- 
gued that the AGN feedback loop can only be established 
if the cooling time is longer than the dynamical time (or 
sound crossing time) of the halo, and showed that allow- 
ing AGN to suppress cooling in hydrostatic haloes resulted 
in a good description of many galaxy formation properties. 
Bow08 extended this by allowing the AGN to expel mate- 
rial from hydrostatic haloes (rather than simply replacing 
the energy radiated). They showed that this model was able 
to match the observed X-ray scaling relations of groups and 
clusters as well as many of the observed properties of galax- 
ies. Typically, the heat input is assumed to be associated 
with low-excitation radio sources (Croton et al. 2006; Best 
et al. 2007), however, the exact heating mechanism is not im- 
portant to the scheme. The crucial distinction is that only 
hydrostatic haloes are affected and that the cold gas disk of 
the galaxy is affected only indirectly because of the reduc- 
tion in the supply of cooling gas (van de Voort et al. (2011) 



discuss these effects in the context of hydrodynamic simula- 
tions). It is important to note that the effectiveness of AGN 
feedback is not a threshold imposed at a fixed mass, but is 
the result of dynamically tracking the relative cooling and 
dynamical times of the halo as it evolves. 

The results of including this type of AGN feedback in 
the SW and MS model are shown in Fig 11. In each model 
we choose the wind speed and mass loading to achieve a 
good match to the abundance of galaxies with stellar mass 
between 10 9 ' 5 and 10 11 /i _1 Mq and have then adjusted the 
Qcooi parameter to achieve a good match to the observed 
stellar mass function. Increasing Q C ooi adjusts the ratio of 
cooling and free-fall times at which the AGN is assumed to 
become effective. Larger values make AGN feedback more ef- 
fective and shift the break in the mass function to lower stel- 
lar mass. With a suitable value for a CO oi the models match 
the observed mass function well. For the SW model, we find 
that (lVind,200)/32oo) ~ (180,8) gives a reasonable match 
to the observed SMF if combined with a C ooi = 0.25. For the 
MS model, we find that (u w jnd,200, P200) = (225, 8) combined 
with a coo i = 0.35 gives a good description of the obser- 
vational data for stellar masses above ~ lO 9 h -1 M . Both 
models over predict the abundance of the lowest mass galax- 
ies, although this problem is reduced in the MS model. To 
emphasise this point, we have superposed pre-liminary data 
from the GAM A survey (Baldry et al., in prep). This pro- 
vides an independent assessment of the mass function for low 
mass galaxies. Although the detailed shape of the mass func- 
tion differs slightly from Li & White (2009), the differences 
are much smaller than the discrepancy between the mod- 
els and the observational data. In contrast the B8W7 model 
keeps a shallow mass function slope down to the faintest 
galaxies plotted (see Fig. 1). The difference in behaviour 
arises from the steeper slope of the stellar mass - halo mass 
relation in B8W7. 

In Fig. 12, we compare the evolution of the mass func- 
tion for the models discussed above. The figure also shows 
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Figure 11. (left panel): Effect of "hot-halo" AGN feedback on the SW model. Both models have (f w i n( j 200; A200) = (180,8), and we 
compare a coo i = 0.0 and 0.25. Increasing « coo i adjusts the ratio of cooling and free-fall times at which AGN are assumed to become 
effective. The increase shifts the break in the mass function to lower stellar mass so that the model matches the break in the mass 
function well, (right panel): We show the effect of AGN feedback in the MS (/3 <x 1/fdiak) model. In this case, we compare models with 
(''wind 200;/^20o) = (225,8) for a coo i = 0.0,0.35. In both panels, black points show observational data. We have supplemented the Bell 
et al. (2003) (circles) and Li & White (2009) (pluses) with preliminary data from the GAMA survey (crosses, Baldry et al., in prep) to 



provide an independent assessment of the observational uncertainties at M, 
low M* is inconsistent with both recent data-sets. 



< 10 h Mq. The up-turn in the model mass functions at 



the B8W7 model, and recent observational data from Drory 
et al. (2005); Bundy, Ellis, & Conselice (2005); Marchesini 
et al. (2009) and Mortlock et al. (2011) (plus, circle, cross 
and triangle respectively). All of the models include AGN 
feedback following the Bow08 scheme. The key issue that we 
wish to test with this plot is whether the models generate 
sufficient massive galaxies at high redshifts, and we focus on 
the brightest galaxies at each epoch. Compared with B8W7 
and the observational data, the new models show a rapid 
decrease in the abundance of the most massive galaxies at 
higher redshift. The discrepancy is worst for the SW model. 

The differences in the behaviour of the models can be 
traced to the differences in wind mass loading in high mass 
haloes (see §3.2.1). The effect arises since the "hot-halo" 
mode of feedback is not a simple cut-off in cooling at high 
halo mass, but explicitly compares the halo cooling time and 
dynamical time, taking into account the halo formation his- 
tory. In practice the effective halo mass threshold increases 
slowly with redshift. Moreover, since none of the models can 
eject material from massive haloes, the efficiency of star for- 
mation is inversely proportional to the mass loading. The 
net result is that massive galaxies appear at higher redshifts 
in the model with the strongest halo mass dependence of 
the mass loading. Since /3 tx in the B8W7 model, this 

model provides the best match to the observational data, 
followed by the MS model (/3 oc v^ lo ). 

We should, however, note that there are considerable 
random uncertainties in determining the stellar masses of 
high redshift galaxies. Applying this convolution will tend 
to smear the model predictions, resulting in a tail of higher 
mass galaxies (see discussion in Marchesini et al. (2009)). 
Thus while this data-set set picks out the B8W7 model, 
careful analysis of the observational errors is required be- 
fore reaching a definitive conclusion. In order to illustrate 
the effect of random uncertainties in the mass determina- 
tion, dotted lines show result of convolving the model with 



a random error of 0.2 dex. This has a pronounced effect on 
the abundance of the most massive galaxies, as a small pop- 
ulation of galaxies that are mistakenly assigned low stellar 
mass can easily overwhelm the true population. The com- 
parison still favours the B8W7 model, but assigning larger 
mass errors would make it difficult to exclude the MS model 
with high confidence. 

Focussing on lower mass galaxies, we see that all of the 
models appear to over-predict the observed normalisation of 
the mass function at z — 1-2 (Marchesini et al. 2009). Al- 
though there is considerable scatter between data-sets, and 
the survey volumes are relatively small, this does appear 
to be a persistent problem, and only the data from Drory 
et al. (2005) are consistent with the evolution seen in the 
models. This discrepancy is also evident if the K-band lu- 
minosity functions are compared directly (eg., Cirasuolo et 
al. (2010)). Pozzetti et al. (2010) suggest that the problem 
lies with the mass dependence of the specific star forma- 
tion rates of the model galaxies, and we will examine this 
in Section 5.1. Our preferred interpretation is that the data 
require a stronger redshift dependence of feedback. We have 
already shown that the pGIMIC model provides a good de- 
scription of the mass function at z > 1, so a promising route 
would be to vary the wind speed parameter between 180 at 
z = and 275 at 2 = 1. Alternatively, the required variation 
in the fraction of the wind escaping would naturally arise 
if we were to choose a criterion for wind escape based on 
the halo mass rather than circular velocity, at least at low 
redshift. It is unclear why this choice should be physically 
motivated, however. Perhaps a better explanation could be 
the greater gas content of high redshift disks, and thus the 
tendency for star formation to occur in more massive star 
forming complexes (Jones et al. 2010; Genzel et al. 2011). 

In summary, introducing a hot halo mode of feedback 
creates a break in the stellar mass function in all three mod- 
els. As a result, all three provide a good match to the ob- 
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Figure 12. The evolution of the SMF compared between three models that match the z = SMF well: B8W7 (red), SW (blue), MS 
(green). All three models include AGN feedback, and the 2 = stellar mass function of B8W7 is repeated as a black dashed line in each 
panel. The effect of observational errors are indicated by the dotted lines which show the effect of convolving the models with a stellar 
mass error of 0.2 dex. For comparison, we show observational data for the high redshift mass function from a variety of sources (see text 
for details) as black symbols. Out to z = 2, all three models result in similar evolution of the stellar mass function. The break in the mass 
function evolves more quickly in the SW and MS models than in the B8W7. However, the deficit of massive galaxies in the MS model 
could be accounted for if the random stellar mass errors are greater than 0.2 dex. All of the models show an excess of M» ~ 10 10 h~ 'Mq 
galaxies at z = 1 — 2 compared to most observational data sets. 



served SMF above a stellar mass of~ 10 9 - 6 A _1 M S . At lower 
stellar masses, the mass function of the SW model rises 
steeply, and is inconsistent with the observational data. This 
trend is less pronounced in the MS model, while the B8W7 
model has a flat SMF to much lower masses. The models pre- 
dict different evolution of the SMF, with B8W7 showing the 
highest abundance of massive galaxies at z = 1 and above. 
All three models predict an abundance of 10 10 /i -1 M galax- 
ies at z = 1 that appears to be at odds with the data, and 
suggest that the effective wind speed should be higher at 
z > 1 than at the present day. 

4.2 The "Starburst" (or "QSO") mode' 

Another channel of AGN feedback, often referred to as the 
"QSO" or "starburst" mode, also has the potential to be 
important because a large fraction of the black hole mass is 
accreted in this way. The GALFORM model assumes that 
black hole growth is triggered when gas is transported to 
the center of a galaxy by disk instabilities or galaxy merg- 



ers. Most of the cold gas fuels a burst of star formation, but 
a small fraction is accreted by the black hole (eg., Springel, 
Di Matteo, & Hernquist (2005)). Since most mergers are gas 
rich, this results in a strong correlation between black hole 
mass and bulge mass very like that observed. We will use 
the term "QSO mode" and "starburst mode" interchange- 
ably, perhaps the term "starburst" should be preferred since 
it makes it clear that this channel only occurs during such 
events. The key distinction is that the "starburst" mode 
acts on the cold gas of the host galaxy, rather than acting 
through the heating of hot gas in the haloes of galaxy groups 
and clusters. In the starburst mode feedback may lead to ex- 
plosive winds that blow cold gas out of the host galaxy. If 
cold gas is removed from the system at sufficiently high spe- 
cific energy it will suffer a long delay before it is able to 
cool once again. This type of feedback has been explored in 
idealised numerical simulations which have shown that the 
energetics of the black hole can plausibly remove the whole 
interstellar medium of the merging galaxies (Springel, Di 
Matteo, & Hernquist 2005; Hopkins et al. 2006) (although 
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Figure 13. (Left Panel) Starting from the SW model 
[(■"wind 200: P200) = (180,8), blue line], we contrast the effect 
of the AGN "hot-halo" mode feedback scheme (implemented by 
setting o coo i = 0.25, green line) with the effect of inducing a 
wind with very high mass loading (,8200 = 16) and wind speed 
(''wind 200 = 1130 kms -1 ) during bursts of star formation (red 
line). We assume that the energy to drive such powerful outflows 
comes from the QSO phase of AGN growth. In order to show the 
effect of more frequent outbursts, we also show an "extreme QSO" 
model in which disks are much more unstable than in B8W7 (cyan 
line), so that black hole fueling occurs more frequently. The main 
effect of the QSO mode is to suppress the formation of galaxies 
near the break of the mass function and to drive the mass function 
towards a power-law shape. The sharp break in the mass function 
is only created if AGN feedback is only effective in hydrostatic 
haloes. 



break in the SMF. We can produce a stronger effect on the 
mass function by increasing the frequency of starbursts. A 
simple way to achieve this is to tighten the disk stability cri- 
terion so that disks more frequently become unstable. The 
effect is illustrated by the "extreme QSO" model in the plot 
(cyan line). The model has been shifted further from the ob- 
served SMF, giving the mass function an almost power-law 
form. 

We can compare these models with Gabor et al. (2011) 
who modify the hydrodynamical models of Op08 to inves- 
tigate the effect of quenching star formation after galaxy 
mergers. They contrast this form of feedback with a model 
in which star formation is suppressed in hot haloes. The im- 
plementation of their schemes are similar to those adopted 
here (although their quasar-mode feedback scheme is trig- 
gered only by mergers, while it is triggered by both merg- 
ers and disk instabilities in our model), and the results are 
qualitatively similar. In particular, the merger model tends 
to have a relatively weak effect on the overall mass function, 
and fails to imprint a characteristic scale on the SMF. 

We have experimented with using other models as a 
starting points. If we start from a model with much weaker 
quiescent feedback, a strong "starburst" mode feedback fails 
to reproduce the shape of the observed SMF, again tending 
to drive the mass function towards a power-law. The "star- 
burst" mode does not have the required effect because it 
does not scale strongly with halo mass (as is the case for the 
"hot-halo" mode). In summary, while the starburst/QSO 
channel might supplement the feedback from supernova dur- 
ing star bursts, it does not provide a scheme for creating a 
break in the stellar mass function. 



higher resolution simulations suggest that the geometry of 
the central outflow may play an important role, Hopkins & 
Elvis (2010)). However, this channel expels only the cold 
material from the system, and does not prevent further ac- 
cretion. As the halo grows, it accretes new satellite galaxies, 
together with their gas, so that (in practice) star formation 
quickly re-establishes itself. 

Fig. 13 illustrates the effect of the QSO mode of expul- 
sion. All of the models we consider reproduce the observed 
correlation between the mass of the black hole and the mass 
of the galaxy bulge. We start from the SW model (with 
parameters (« w ind,20o, P200) = (180,8), blue line). In this 
model, black holes grow strongly as a result of galaxy merg- 
ers and disk instabilities (see Bow06), but this results in no 
effective feedback. The default model assumes that all the 
energy generated in black hole events is radiated without 
doing significant mechanical work. In order to explore what 
would happen if this radiation coupled effectively to the sur- 
rounding gas (or if the quasar accretion disk produced a high 
speed wind), we implement a "QSO mode" of feedback by 
using much stronger feedback during starbursts (compared 
to quiescent star formation events). We illustrate the effect 
by using (« W ind,200, P200) = (1130, 16) during bursts (the re- 
sults are similar for other parameter choices) so that the 
energy of the wind during the burst is 80 times larger than 
that during quiescent star formation. The figure shows that 
even winds of this strength have a modest effect on the mass 
function. Furthermore, their effect is to suppress the abun- 
dance of M* galaxies rather than to create an exponential 



5 FURTHER CONSIDERATIONS 
5.1 The Star Forming Sequence 

We have compared the different schemes on the basis of 
the stellar mass function and the central stellar mass. In 
this section, we compare the models with the star formation 
rates of the central galaxies. We focus on B8W7 and the best 
fitting models with constant wind and /3 oc 1/udisk feedback 
scalings, SW+AGN [(w wtod ,2oo, #200, <w) = (180,8,0.25)] 
and MS+AGN [(iWd.aoo, P200, a coo i) = (225,8,0.35)]. The 
SMFs of these models are shown in Fig. 11. 

Fig. 14 shows the logarithm of the z = specific star 
formation rate (SSFR) as a function of stellar mass. We 
focus on the properties of star forming galaxies (which we 
define as having SSFR > 0.01 Gyr -1 ). In the model we re- 
strict attention to central galaxies to avoid uncertainties in 
the treatment of satellite galaxies. In all the models, there is 
a clear sequence of star forming galaxies that can be cleanly 
compared to the observed star forming sequence. There is 
also a large population of galaxies which are not seen on this 
diagram because their star formation rates are extremely 
low. These are satellite galaxies, or central galaxies in haloes 
with effective AGN feedback. 

We compare the theoretical models to observational 
data from Brinchmann et al. (2004, updated to DR7). The 
observational data is shown as black contour lines in the 
figure. The observed relation is almost flat (with low mass 
galaxies having slightly higher SSFR) up to stellar masses 
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Figure 14. Comparison of the specific star formation rates of 
central galaxies for three models which provide encouraging fits 
to the SMF. In seperate panels, we show the baseline model, 
B8W7, SW+AGN [(« win d,200 , £200 , <* cool ) = (180,8,0.25)] and 
the MS+AGN model [K ind , 20 o, P200, <*cool) = (250,8,0.35)]. 
The star formation sequence for local galaxies is shown as black 
contours (see text for details). A random sample of model galax- 
ies are shown as blue points, with the median relation and 10% 
and 90% percentiles shown as a red line with error bars. The 
B8W7 model has a constant specific star formation rate regard- 
less of stellar mass. In contrast, the SW and MS models produce 
relations that are strongly dependent on system mass and are 
incompatible with the observational data: instead of a near con- 
stant specific star formation rate, the SW and MS model have 
specific star formation rates that are lower in small galaxies. 



of 10 11 /i _1 M0 (above which star formation is suppressed 
by AGN feedback). Given the age of the universe, the loca- 
tion of the sequence at SSFR ~ 10~ ' 8 Gyr~* implies that 
the specific star formation rates of all galaxies have been 
steady (or slowly rising) over the history of the universe. 
Note that the uncertainties due to dust obscuration would 
tend to increase the star formation rates of the most mas- 
sive galaxies. However, the data presented already include 
an extinction correction based on the Balmer decrement, 
and agree well with specific star formation rates based on 
SED fitting (McGee et al. 2011). 

In each panel, a random sample of model galaxies are 
shown as blue points, while the red line with error bars shows 
the median specific star formation rate and the 10 th and 90 th 
percentiles of the distribution. We include only central star 
forming galaxies (with SSFR > 0.01 Gyr -1 ) in this calcula- 
tion, but very similar results are obtained if we include star 
forming satellite galaxies as well. The observed relationship 
is reproduced fairly well by B8W7. The specific star forma- 
tion rate in the model is flat over a wide range in stellar mass, 
from below 10 7 h~ 1 M e to almost ~ 10 11 h~ 1 M Q (where the 
relation dips as the supply of fuel for star formation is sup- 
pressed by the hot-halo feedback). 

In contrast, the SW and MS models predict a rela- 
tion with a noticeable decline in SSFR towards lower stel- 
lar masses. Although the declining trend is less abrupt in 
the MS model (bottom panel) than in SW (middle panel), 
both relations are clearly inconsistent with the observational 
data. We infer from this figure that a dramatic change of 
feedback efficiency cannot be responsible for the flattening 
of the SMF. This problem is also seen in the hydrodynamical 
models of GIMIC and Op08 (see Dave, Oppenheimer, & Fin- 
lator (2011)), although the limited mass resolution of those 
simulations limits the comparison to galaxies more massive 
than ~ 10 9 /i _1 M0, and consequently the mass dependence 
is not so clearly evident. 

It is interesting to understand the origin of the dip. In 
the SW and MS models, the flat region of the SMF is created 
by a transition between the two feedback regimes: for low 
halo masses, feedback is extremely effective at expelling gas 
from the halo and the star formation rate is strongly sup- 
pressed. At higher masses, however, the wind velocity is no 
longer sufficient to escape the halo and the cold gas mass and 
star formation rate increase. However, because the division 
between the two regimes occurs at a fixed escape velocity, we 
expect that the halo mass of the transition evolves rapidly 
with redshift, Mhaio, kink ~ (1 + z)~ 3 ^ 2 . Allowing for the de- 
pendence of stellar mass on halo mass, M* tx Mhaio 2 (eg. 
Fig. 2, the relation evolves slowly with redshift) we expect 
the stellar mass at which the transition occurs to evolve as 
Af„,kink ~ (1 + z)~ 3 . Thus the transition mass evolves more 
quickly than the mass of an individual galaxy. Thus the 
transition mass is much smaller at high redshift. Over time, 
an individual galaxy makes a transition from the regime 
in which ejection is in-effective to the one in which it is. 
Consequently, galaxies somewhat below the transition mass 
at z = have low current star formation rates compared 
to their past average. At the very lowest stellar masses, the 
SSFR in the SW and MS models begins to increase. Galaxies 
that lie well below the kink in the M* — Mhaio relation have 
experienced similar feedback during their formation history 
and the rise is thus to be expected. 
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Although the B8W7 model provides the best descrip- 
tion of the observational data, it does not reproduce the 
weak trend for the SSFR to increase as the stellar mass de- 
creases (SSFR oc A/*" 22 ) that is seen in the data. The 
strength of this trend is controversial, but does not appear 
to be an observational selection effect. It is seen regardless 
of the star formation diagnostic that is applied, and is ap- 
parent at higher redshifts as well as locally (but this de- 
pends critically on the definition of the sample — for a re- 
cent overview, see Karim et al. (2011)). In §4.1, we noted 
that the surprisingly rapid evolution of the normalisation 
of the observed mass function suggested that the effective 
wind speed should scale with redshift. This change would 
also have implications for the SSFR, since the present-day 
star formation rate would increase relative to the past aver- 
age. Varying the escape speed rather than the mass loading 
could create a tilt in the SSFR — M* relation since the effect 
will be strongest around the kink in the M* — M^aio relation 
but result in little change in the star formation histories of 
the most massive star forming galaxies. 

In summary, the specific star formation rate of galaxies 
provides an important additional discriminant of the models. 
The B8W7 model comes closest to matching the observed 
data, with the characteristic specific star formation rate that 
is almost independent of stellar mass. In contrast the SW 
and MS models show specific star formation rates that de- 
cline with decreasing stellar mass, while the observational 
data show an slightly increasing trend. Further exploration 
of feedback schemes that scale systematically with redshift is 
required to identify a model which produces a better match 
to the observational data. 



5.2 The Star Formation History of the Universe 

We have seen that all the models reproduce the observed 
build-up of stellar mass reasonably well, another way to 
tackle this question is to investigate the star formation rates 
of galaxies directly. In Fig. 15 we show the evolution of the 
cosmic star formation rate. For the models this is calculated 
by integrating the contribution of galaxies down to stellar 
masses of 10 s /i _1 Mq. As a result of including feedback from 
AGN in order to suppress the formation of the most massive 
galaxies, the behaviour of all three models is broadly similar 
in this plot. Above z = 2, the integrated star formation rates 
of the models are very similar. However, the strength of the 
decline in the star formation rate differs between the mod- 
els, being much stronger for the B8W7 model than for the 
SW scheme. As expected, the MS scheme lies in between. 

We compare the models with the compilation of ob- 
servational data from Hopkins (2004) and the more recent 
data of Karim et al. (2011) (radio, red circles), Rodighiero et 
al. (2010) (Mid-IR, blue triangles) and Cucciati et al. (2011) 
(rest-frame UV). The radio and mid-IR based measurements 
have the advantage of that the obscured star formation is ac- 
counted for (see Karim et al. (2011) for further discussion). 
None of the models match the trends in the observational 
data perfectly. Given the scatter in the observational data 
sets there is little reason to choose between the MS and 
B8W7 models on the basis of this plot. However, the uncer- 
tainties in the observational data arise in large part because 
of the large extrapolation required to correct the observed 
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Figure 15. The evolution of the cosmic star formation rate den- 
sity as a function of redshift. Solid lines show the total star for- 
mation rate density. Line colours distinguish different feedback 
models. All three models include hot-halo AGN feedback that 
suppresses the formation of stars in high mass haloes. While the 
models all show similar star formation rates above z = 2, the 
decline in the star formation rate between 2 = 1 and the present 
day differs markedly. The solid lines should be compared with 
the observational data, shown as black crosses (from the com- 
pilation of Hopkins (2004)), red circles (from the stacked VLA 
analysis of Karim et al. (2011)), blue triangles from Rodighiero 
et al. (2010) and black triangles from Cucciati et al. (2011). The 
dashed and dotted lines show the contribution from galaxies more 
massive that 10 9 ' 5 and 10 10 ' 5 h~ 1 M.Q respectively. Despite overall 
the similaraity of the star formation histories in the three models, 
the way the star formation rate is divided between stellar masses 
varies greately. 



star formation rate density for galaxies that are too faint to 
be directly detected. 

An important question is therefore to examine how the 
star formation rate density depends on the stellar mass of 
galaxies. The mass dependence of the cosmic star formation 
rate in the models is illustrated by dashed and dotted lines in 
Fig 15. The line-styles show the integrated contribution from 
galaxies more massive than 10 9 ' 5 and 10 10 ' 5 /i _1 Mq respec- 
tively. Although the total star formation rates of the models 
are similar, the way the star formation is distributed be- 
tween galaxy masses is different, with the SW and MS mod- 
els showing a shift from star formation that is dominated by 
low mass (< 10 9 ' 5 /["'Mq) galaxies at high redshift (note 
the large difference between the solid and dashed curves) to 
a dominance by lower mass galaxies at low redshift. In con- 
trast, in the B8W7 model massive galaxies make a similar 
contribution to the total star formation rate density at all 
redshifts (ie., there is a constant offset between the solid, 
dashed and dotted lines). Thus, although all three models 
show a similar total star formation rate density at z > 2, 
the contributions of different mass galaxies is very different 
in the models. 

Given the different contributions of high and low mass 
galaxies in the models, it is very dangerous to quantitatively 
compare to observational data for the total star formation 
rate density since the observations are usually based on ex- 
trapolation of the properties of high mass galaxies. There- 
fore in Fig. 16 we show the mass dependence of the cos- 
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Figure 16. The panels in this figure show the evolution of the 
cosmic star formation rate density as a function of stellar mass. 
The coloured lines differentiate different feedback schemes, while 
the panels compare to data at z = 0.0, 1.0 and 2.2. Solid lines 
show the distribution obtained from the models, dashed lines il- 
lustrate the effect of a random error of 0.2 dex in the stellar mass 
estimates. Observational data collected by Gilbank et al. (2011) 
and Karim et al. (2011) are shown as black symbols. This figure 
gives a more accurate comparison of the observed star formation 
rate density with that predicted by the models (see text for dis- 
cussion). 



mic star formation rate, separating the redshift dependence 
into separate panels. The different models are distinguished 
by the coloured lines. The solid lines show the distribution 
measured for the model galaxies, and the dashed lines illus- 
trate the effect of a 0.2 dex error in stellar mass assignment. 
Using this plot, it is then possible to directly compare to 
observational measurements based on mass-complete sam- 
ples (Gilbank et al. 2011; Karim et al. 2011). Historically, 
such plots have been use to infer that the star formation 
density is dominated by large galaxies (at high redshifts), 
and by lower mass galaxies at low redshift. In practice, the 
more complete and sufficiently deep data sets show that the 
increase in the star formation rate density with redshift is 
similar for all galaxy masses, however, the limitations of the 
observed data-sets, particularly at high redshift are evident. 

Comparison of the models and observational data in 
this logarithmic plot allows the contribution from the tails of 
the mass distribution to be clearly seen. At low redshift, star 
formation in the SW and MS models is more concentrated 
towards high mass galaxies than suggested by the data. In 
contrast, star formation in the B8W7 model is somewhat 
too flat. Also, while a significant contribution to the star 
formation rate comes from low mass galaxies (as suggested 
by the observations), the model fails to generate sufficient 
massive star forming galaxies. This occurs because the AGN 
feedback parameters need to match the observed SMF set 
an effective halo mass threshold that is lower in B8W7 than 
the SW and MS models. This is driven by the larger scatter 
in galaxy mass at fixed halo mass in B8W7. Thus the model 
contains galaxies of the stellar mass required to match the 
observed SMF, but star formation has been suppressed in 
too many of them. A stellar mass error of 0.2 dex improves 
the match (as shown by the dashed line), but this probably 
over-estimates the uncertainty in the local data. The differ- 
ence in behaviour at z = tallies with the differences in the 
models' SSFR-M* relations seen in Fig. 14. 

The behaviour of the models is similar at z = 1: B8W7 
results in a distribution that is flatter than the data, while 
the MS and SW models are slightly too peaked. The or- 
der of the models at the massive end is reverse, with the 
B8W7 model suggesting the greater population of massive 
star forming galaxies. The differences are small, however, 
and likely to be masked by the uncertainties in the stellar 
mass determination (see dashed lines). At z = 2, a deficit of 
massive star forming galaxies is apparent for the SW model, 
and the rapid drop at high masses is not reconciled with the 
observational data if an error 0.2 dex is assigned to the stel- 
lar mass. This deficit is to be expected from our comparison 
with the SMF of this model: the SW model lacks sufficient 
high mass galaxies. 

In summary, while all models result in similar predic- 
tions for the evolution of the integrated cosmic star forma- 
tion rate, the mass dependence of the cosmic star forma- 
tion rate exposes important differences between the mod- 
els. Although none of the trends are conclusive, the com- 
parison highlights different issues with the models. The SW 
and MS models tend to under-represent star formation from 
low mass galaxies at low redshift. At high redshift, the SW 
model under-represents the contribution to star formation 
from high mass galaxies. In contrast, the B8W7 tends to 
under-represent the contribution from high mass galaxies at 
lower redshift. As the consistency and completeness of the 
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Figure 17. The fraction of baryons (relative to Mhaio^W^o) i n 
different phases as a function of halo mass. The lines are cumula- 
tive showing, blue: stars, green: adding cold gas, red: adding hot 
and reheated gas. The three panels show the B8W7, SW and MS 
models (all including the "hot-halo" mode of AGN feedback). The 
models include the effect of hot-halo AGN feedback, as shown in 
Fig. 11. The shaded region indicates the total stellar mass con- 
tent of haloes measured in observational data-sets (see text for 
details). 



observational data improves this approach has great poten- 
tial for discriminating between feedback models. 



5.3 Where are the baryons? 

As we have seen, in the successful models, only a small frac- 
tion of baryons are locked up into stars. It is of interest to 
examine the phase of the remaining baryons in the model, 
and this is illustrated in Fig. 17. We focus on the same three 



models that provided a good description of the SMF above 
M» ~ IO^^Mq: B8W7, SW (tWd.200, £200) = (200,4) 
and MS (lWd,200, P200) = (200,16). The SMFs of these 
models are shown in Fig. 11. 

The distribution of baryonic matter between phases is 
plotted as a function of halo mass. We include all baryons 
associated with the halo by the model, and define the baryon 
fraction as the mass of baryons, Mb divided by the cosmic 
average (ie., ft, = Mj,/f2;,Mhaio). In the absence of feedback, 
the baryons collapse with the dark matter, and all haloes 
would be baryonically 'closed' with fractions of order unity 
(Crain et al. 2007). The lines are cumulative so that the 
mass in stars is shown as a blue line. The green line adds 
the contribution of cold gas, and the red line includes the 
contribution from diffuse halo gas (ie., 'hot' plus 'reheated' 
components in the terminology of Bow06) This gives the 
total baryon mass associated with the halo. Note that, be- 
cause of the expulsion feedback considered in these models, 
the total mass of baryons in the halo need not reach the 
cosmic value. 

As we should expect, given that the models reproduce 
the M* — Mhaio relations deduced from sub-halo abundance 
matching, all the models have a similar peak in the fraction 
of baryons locked into stars around Mhaio ~ 10 11 ' 5 h~ Mq 
(solid lines). In the B8W7 model, this declines rapidly to 
lower masses, and declines to an almost constant value 
(~ 0.08) in higher mass systems. In the SW model, the ex- 
cess abundance of small galaxies means that the decline to 
lower mass haloes is less steep, reaching a constant value of 
0.05 in haloes smaller than 10 11 h~ 1 M Q . The MS model is 
intermediate between the two, as we should expect. The am- 
plitude of the peak is smaller in the B8W7 case than in the 
other models. These differences are largely due to the contri- 
bution from very faint satellite galaxies (M» < 10 9 ^'Mq) 
in the SW and MS models, although it is also noticeable that 
the peak in stellar abundance is broader in these systems. 
This is consistent with all of the models giving a similar 
match to the mass function because of the interaction be- 
tween the halo abundance and the scatter in the M* -Mhaio 
relation. 

The most striking difference between the plots lies 
in the behaviour of the total baryon contributions below 
Mhaio ~ 10 11 ' 5 ft _1 M . In the B8W7 model (in this regime), 
all of the baryons reside in the host halo, and the diffuse gas 
makes up the majority of the baryons. It is questionable 
whether such a large halo component (which we refer to 
as the 'circum-galactic medium', CGM) is compatible with 
observational data. Even though recent observations from 
the Cosmic Origins Spectrograph (COS) on HST suggest 
that there is a significant CGM (Tumlinson et al. 2011), 
it is likely to account for "only" a smilar fraction to the 
mass in the galaxy's stars (ie., ~10% of see Prochaska 
et al. (2011) for discussion). In contrast, the SW and MS 
models are able to eject much of the material from the halo 
so that low mass galaxies contain little diffuse gas in their 
haloes, and most of the baryons are locked into cold gas. As 
cold gas is ejected from these galaxies it escapes from the 
halo. Nevertheless, the feedback scheme in these models has 
been tuned to match the observed stellar mass function, so 
that the outflowing wind stalls as the halo mass approaches 
10 12 /i _1 Mg, resulting in a peak in the abundance of diffuse 
gas that matches the peak in the stellar fraction. At masses 
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above ~ 10 12 /i _1 M , material is driven out of the halo by 
AGN feedback in all three models, resulting in a second dip 
in the abundance of diffuse gas. As discussed in Bow08, this 
results in a good match to the X-ray luminosities of groups 
and clusters. The rapidly rising gas fractions accounting for 
the tilt in X-ray scaling relations compared to the 'self sim- 
ilar' predictions. 

In order to compare the models with observational data, 
we show measurements of the stellar mass fraction of haloes 
from Leauthaud et al. (2011) as a shaded region in Fig. 17. 
This analysis is based on combining a sub-halo occupa- 
tion distribution model and integrated stellar mass measure- 
ments for a sample of COSMOS X-ray detected groups, for 
which halo masses have been determined using a weak lens- 
ing analysis. In contrast to some previous results (eg., Gio- 
dini et al. (2009); Gonzalez, Zaritsky, & Zabludoff (2007)), 
their careful analysis shows that the stellar mass content of 
high mass haloes is low and comparable to the mean stellar 
mass fraction of the universe as a whole (as would be ex- 
pected on quite general grounds (Balogh et al. 2008)). The 
shaded region indicates the plausible systematic uncertain- 
ties in the analysis. The survey includes only directly de- 
tected galaxies, and could under-estimate the stellar mass 
content of some haloes with significant intra-cluster light. 
However, the correction for intra-cluster light is likely to 
be small and certainly less than 50% (Zibetti et al. 2005; 
McGee & Balogh 2010). As Leauthaud et al. (2011) discuss 
the differences from previous work arises from the treat- 
ment of stellar populations and satellite galaxies, not from 
the contribution of the intra-cluster light. 

All of the models compare with this observational data 
reasonably well. The SW and MS models fit particularly 
well, with a peak in the stellar mass fraction close to 
~ 1O 12 /i _1 M . At higher masses, the observations show a 
shallow decline, which is also seen in the models. Differ- 
ences between the model and the data could certainly be 
accounted for by the missing intra-cluster light distribution. 
The B8W7 model show a somewhat sharper decline than is 
evident in the data, together with a peak in the stellar mass 
fraction that is shifted to lower masses. However, this part of 
the observational region is dominated by the scatter in the 
M*-M ha i relation: all of the models presented reproduce 
the observed SMF, as seen in Fig. 11. 

In summary, the different feedback schemes result in 
significant differences in the distribution of mass between 
the phases of baryonic matter. In the B8W7 model, haloes 
below ~ 10 11 ' 5 /i _1 M s retain the cosmic baryon fraction, 
mostly in the form of a diffuse gas halo. In the SW and 
MS models, most baryons have been ejected from low mass 
haloes. In higher mass haloes, all three models suggest that 
the stellar fraction should quickly converge to the average 
value. The shape of the models fits well with observational 
measurements. 

5.4 Metal Enrichment of the Inter-Galactic 
Medium 

Although the models present have all been developed in or- 
der to explain the galaxy population of the universe and the 
X-ray scaling relations of galaxy clusters, the metal enrich- 
ment of the intergalactic medium could offer an important 
discriminator between the models. Indeed a major motiva- 



tion for considering powerful winds is to explain the wide- 
spread metal enrichment of the Universe. Recent observa- 
tions have shown that metals are widely distributed in the 
inter-galactic medium. For example, Prochaska et al. (2011) 
show that the low redshift metal absorption line cross- 
section is compatible with 0.3 L* galaxies being surrounded 
by a metal enriched halo extending out to ~ 300 kpc (ie., 
well beyond the virial radius of the galaxy). Quantitative 
comparison with these data is fraught with difficulties, how- 
ever, since the observations must be corrected for the ion- 
isation state of the absorbing clouds. This issue is further 
complicated by the multiphase nature of the absorbers. The 
inherent correlation of bright and faint galaxies must also be 
taken into account, particularly if smaller galaxies dominate 
the enrichment. In addition, the phenomenological mod- 
els include a number of simplifying assumptions that make 
quantitative comparison with the data difficult. In particu- 
lar, the models we have presented here assume that enrich- 
ment and metal re-cycling is instantaneous and that metals 
are distributed uniformly through the galaxy's halo. These 
limitations must be carefully taken into account in a detailed 
comparison, and we will not attempt this here. Nevertheless, 
it is instructive to briefly examine the differences between 
one model and another. 

All of the models result in significant ejection of metals 
from the galaxy disk, and there is little to distinguish be- 
tween the models on the basis of the mass of metals ejected 
from the disk. Rather, the important discriminator is the 
distribution of the metals. In the B8W7 models the metals 
remain trapped in the dark matter potential of the galaxy, In 
the SW and MS models, metals ejected from small galaxies 
will be more widely distributed. At face value, the observed 
wide-spread distribution of metals would favour the SW or 
MS models. 

In practice, however, the situation is more complex. At 
low redshift, in all of these models, most of the intergalactic 
metals are associated with material that has been ejected 
from galaxy groups by the action of AGN feedback. The sit- 
uation at low redshift is illustrated in Fig. 18. The panels on 
the left hand-side show the distribution of the baryon mass 
between haloes of different mass, such that the integral un- 
der the curve gives the total baryonic mass of the universe in 
haloes more massive than 10 10 /i -1 Mq. The coloured lines il- 
lustrate the division of the baryons between different phases 
for the different models. The blue, green and red lines show 
the cumulative contribution from star, cold gas and diffuse 
gas, as previously discussed. The black line shows the contri- 
bution to the baryon mass in the absence of expulsion. It is 
noticeable that this is widely distributed, with all haloes less 
massive than 10 14 h~ 1 M & making a similar contribution to 
the total mass content. Because the distribution of mass be- 
tween different haloes is rather flat, the coloured lines bear 
close resemblance to those in Fig. 17, and attention can be 
focused on the difference between the red and black lines. 
Since the scale is linear, this is the mass contribution from 
material that has been ejected from haloes. In the B8W7 
model, this contribution only comes from the material ex- 
pelled by AGN in haloes more massive than 10 12 /i~ 1 M Q , 
while in the SW and MS models, expulsion feedback means 
that there is a second important contribution from haloes 
less massive than 10 11 ' 5 . The integral of the difference be- 
tween the red and black curves measures the total mass frac- 
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Figure 18. (left-hand panels): The contribution of baryons in different phases to the total baryon content of the universe at z = 0, 
plotted as a function of halo mass. The lines are cumulative showing, blue: stars, green: adding cold gas, red: adding hot and reheated 
gas, and can be compared to the mass fraction (per halo mass bin) in different phases shown in Fig. 17. The solid black line shows the 
total contribution associated with haloes of a given mass in the absence of any mass loss. The difference between the red and black lines 
indicates the contribution to the total intergalactic baryon mass that has been expelled from haloes, (right-hand panels): the contribution 
of haloes and phases to the total metal mass content of the universe at z = 0. The black line indicates the metal mass that would be 
associated with haloes in the absence of mass ejection, so that the difference between the black and red lines indicates the contribution 
to intergalactic metals. The figure shows that the contribution to intergalactic metals is dominated by the mass expelled from galaxy 
groups by AGN feedback in all models. This dominates over the metals that are expelled from the haloes of low mass galaxies. 



tion ejected from haloes. We find that the fractions are 52%, 
61% and 66%, for the B8W7, SW and MS models respec- 
tively (for haloes more massive than IO^/i-^q). 

The panels on the right show the equivalent plot for 
the contribution to the metal content of the Universe. The 
line colouring follows that used in the left-hand panels. The 
black line shows the contribution from stars that are now 
in haloes whose masses are given on the horizontal axis. Of 



course, when the metals were produced, the stars may have 
been in lower mass haloes. The metals locked into stars, cold 
gas and diffuse gas are shown by the coloured lines, and the 
area under the curve gives the total contribution. In con- 
trast to the left hand panels, the metal contributions are 
much more strongly weighted to haloes more massive than 
10 12 h~ 1 MQ, which reflects the low efficiency of star forma- 
tion in smaller haloes, and thus the low rate of metal produc- 
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tion. As a result, the contribution to the global metal bud- 
gets from low mass haloes is strongly down-weighted com- 
pared to the contribution from galaxy-group haloes (driven 
by AGN feedback, see below). The fraction of metals that 
have been lost from haloes is given by the difference be- 
tween the red and black lines. To be clear, the difference 
between the lines corresponds to the fraction of metals that 
are absent from haloes of the mass given on the x-axis at 
z — 0. Some of these metals may have been ejected from 
projenitor haloes of smaller mass. We find that the mod- 
els have total intergalactic metal fractions of 42%, 45% and 
46% for the B8W7, SW and MS models respectively. In the 
B8W7, the black and red curves are superposed below a halo 
mass of 10 11 ' 5 /i _1 Mq, but the results are similar for all the 
models, with by far the largest contribution coming from 
the group-scale haloes. Above 10 12 /i _1 M© the behaviour of 
the models is similar since all the models share a common 
AGN feedback scheme. In the absense of AGN feedback very 
few metals are ejected from such haloes. Focussing on lower 
mass haloes, greater differences are apparent. In particular, 
the SW and MS models have an extended tail making a use- 
ful contribution to the intergalactic metals. This tail is not 
present in the B8W7 model, both because the material re- 
mains trapped in the halo potential and because the SMF 
is much flatter in the B8W7 model than in either the SW or 
MS models. Nevertheless, low mass haloes make a relatively 
little contribution to the total metal content of the universe 
in all the models. 

At higher redshifts, the fraction of intergalactic metals 
drops rapidly in all the models, reflecting the lower abun- 
dance of high mass haloes with effective "hot mode" feed- 
back. At these redshifts, a larger fraction of metals come 
from dwarf galaxies in low mass halos and this is reflected 
in greater differences in the inter-galactic metal fraction be- 
tween the models. In Fig. 19, we show the evolution of the 
fraction of metals that have been ejected from galaxies. To 
be specific, we show the fraction of all metals produced at 
higher redshifts, associating ejected metals with the halo in 
which the stars that produced the metals are located at the 
specified redshift. Metals produced by high redshift galaxies 
may be recaptured by larger haloes at low redshift, just as 
the original stars of the early galaxy may have been merg- 
ered into a larger system. In principle, metals that are miss- 
ing from haloes of mass greater than 10 12 /i _1 M© may have 
been ejected at earlier time from a progenitor halo of much 
lower mass; however, by re-running the models without AGN 
feedback, we have shown that this process is negligible. In 
haloes more massive than 1O 12 /i -1 M0, the loss of metals 
from the halo is entirely due to the effect of AGN feedback. 

The black line shows the total fraction of metals that are 
outside galaxy haloes, while the coloured lines show how this 
fraction depends on halo mass: the blue (red) line shows the 
material that has been lost from haloes less massive than 
lO 11 h" 1 M (1O 12 /i _1 M ) so that the contribution from 
material ejected from small (large) galaxies can be seen. 
The blue line is not visible for the B8W7 model because 
material is not expelled from galaxy haloes (the small con- 
tribution from the green line is due to hot-halo mode feed- 
back in Mh a i ~ 10 12 /i _1 Mq systems). The area above the 
red line shows the contribution from metals ejected by the 
"hot halo" mode of feedback. Although this fraction declines 
quickly with increasing redshift it dominates the galaxy con- 
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Figure 19. The fraction of metals that have been expelled from 
haloes as a function of redshift for the three different models. 
The black line shows the total fraction of ejected metals, while 
the red and blue lines show the fraction ejected from haloes less 
massive than 10 12 (["'Mq and 10 11 h~ 1 M.Q respectively. At low 
redshift the inter-galactic metal content is dominated by ejection 
from galaxy group and cluster haloes by the "hot halo" mode 
of feedback. At high redshift there are few such massive haloes 
this ejection mechanism is unimportant. As a result there are few 
intergalactic metals in the B8W7 model at z > 3. In the SW and 
MS models, the fraction of metals ejected from haloes is much 
higher at this redshift, primarily as the result of ejection from 
dwarf galaxies. 



tribution out to z ~ 2, even in the MS and SW models. At 
z > 3, the IGM metal distribution reflects the ejection of 
metals from dwarf galaxy haloes in these models, and the 
observed widespread enrichment of the high redshift IGM 
(Aguirre et al. 2008; Frank et al. 2010; Steidel et al. 2011) 
seems to favour the expulsion models (SW and MS) over 
B8W7 (Bertone, Stoehr, & White 2005; Booth et al. 2010; 
Cen & Chisari 2011). 

Although these figures make for an interesting compar- 
ison between models, quantitative comparison to the ob- 
served abundance of intergalactic metals is extremely diffi- 
cult. In particular, the current model does not specify how 
far outside the group haloes the metal distribution will ex- 
tend, and does not begin to address its temperature and ion- 
isation state. This issue can only be satisfactorily addressed 
if we account for the non-spherical nature of halo accretion 
and outflow. Equally, the extended diffuse gas distribution 
within haloes will also make a contribution to the metals 
seen in random sight lines. It would therefore be premature 
to rule out any of the models without careful consideration 
of the effects of galaxy clustering. This task is beyond the 
scope of the present paper but is clearly an important avenue 
for future work. 

In summary, all three models that we have considered 
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eject the majority of metals from the galaxy disk. The met- 
als locked into stars and cold gas represent a relatively minor 
component of the universe. In the expulsion models, SW and 
MS, the majority of metals completely escape from haloes 
smaller than 10 11 ' 2 h~ 1 M.Q. Thus these metals are widely 
distributed across intergalactic space. In contrast, the B8W7 
model retains metals within lower mass haloes, creating a 
circumgalactic medium. However, the dominant contribu- 
tion to intergalactic metals in all three models comes from 
the material ejected from group-scale haloes by the action 
of AGN feedback. While current observational data suggests 
that the metal haloes of galaxies extend well beyond their 
virial radii (and thus that the haloes are more extended 
than the B8W7 model would suggest for individual galax- 
ies), a more detailed study of the effects of clustering (and 
the observed clustering of metal lines) is required before the 
models can be distinguished on this basis. 



6 CONCLUSIONS 

'Feedback' is a fundamental component of galaxy formation 
models, allowing us to understand the marked differences 
between the dark matter halo mass function and the galaxy 
stellar mass function (SMF). In order to explain the obser- 
vational data, galaxy winds must strongly suppress the for- 
mation of low-mass galaxies. Phenomenological (or "semi- 
analytic") models, such as GALFORM, have been shown to 
provide a good description of the galaxy mass function to- 
gether with its evolution, the specific star formation rates of 
galaxies and the contribution of galaxies of different masses 
to the cosmic star formation rate density. In the standard 
GALFORM model, these successes are achieved by adopt- 
ing a feedback parameterisation that varies strongly with 
system mass, such that the speed of the outflow tracks the 
mass of the halo while the mass loading of the wind decreases 
strongly with halo mass. In contrast most hydrodynamical 
models, including the GIMIC simulations, adopt a wind that 
is independent of the system mass. 

In this paper, we have modified the semi-analytic code 
to implement feedback schemes similar to those usually 
adopted in hydrodynamic models. Our scheme allows the 
wind speed and mass loading to be adjusted and the con- 
sequent effect on the stellar mass function of galaxies, its 
evolution and other galaxy properties to be determined. We 
focus on three particular models. Two models have fixed 
wind parameters model: the pGIMIC model has high wind 
speed and modest mass loading similar to the GIMIC hy- 
drodynamical simulations (Crain et al. 2009); the SW model 
has a slower wind speed (at the same mass loading). We also 
consider a model (MS) in which the wind mass loading scales 
with the inverse of the disk circular velocity. This is similar 
to the momentum scaling models of Oppenheimer & Dave 
(2008), although we set the wind speed by keeping the total 
energy of the wind (rather than its total momentum) fixed 
across halo masses. Initially, we also considered a model with 
a l/w disk dependence of mass loading (ES), but found this 
could not reproduce the observed stellar mass function. We 
contrast these models with a model based on Bow06 and 
Bow08, but adapted to the WMAP7 cosmology used in this 
paper (B8W7). For each of these models we consider the 



effect of including a "hot-halo" mode of AGN feedback fol- 
lowing the gas expulsion scheme of Bow08. 
The main results are as follows: 

• We find that the phenomenological description of feed- 
back that we use in our code reproduces many of the trends 
seen in hydrodynamic simulations. This give us confidence 
that the phenomenological approach captures the key physi- 
cal effects of galactic winds well and allows us to explore the 
parameter space of winds properties quickly and efficiently. 
Comparison of the phenomenological model and hydrody- 
namic simulation also allows us to better understand the 
key physical processes in galaxy formation. A model which 
uses a high wind speed and moderate mass loading (and does 
not include AGN feedback), reproduces the form of the SMF 
seen in the GIMIC numerical simulations well. Although 
the strong feedback suppresses the stellar masses of galaxies 
strongly, the effect is most prominent around 10 10 ' 5 ft _1 MQ, 
so that the resulting SMF has a dip in the abundance of 
galaxies at these masses. At lower masses, the SMF rises 
steeply since star formation is too weakly suppressed, while, 
at higher masses the wind is not sufficiently energetic to es- 
cape the halo. Using the phenomenological model, we are 
able to explore how the dip in the mass function depends on 
the assumed wind speed and mass loading. 

• The observed mass function is better matched by adopt- 
ing a relatively slow wind, as seen in the hydrodynamic 
simulations of Op08. The SW model, with a wind speed of 
180 km s -1 and a wind mass loading of 8, results in a good 
match to the observed galaxy mass function over the range 
10 9 - 10 11 /i _1 M . Introducing the AGN scheme developed 
in Bow08 suppresses the formation of galaxies from cooling 
in hot gas haloes, resulting in a good match to the high-mass 
turn-over of the mass function. The combined (AGN plus 
supernova-driven wind) model provides a good description 
of galaxy abundance down to 10 9 /i~ 1 M Q (but rises much 
more steeply than the observational data at lower masses). 
Although this is encouraging, we find that the model per- 
forms poorly compared to B8W7 in several other aspects. 
In particular, when we examine the specific star formation 
rates in this model, we see that the fixed feedback scheme 
results in a kink in the specific-star formation rates of galax- 
ies at around a stellar mass of 10 9 ' 5 /i- 1 M Q , such that low 
mass galaxies have a factor ~ 5 lower star formation rates 
than their massive counter-parts. This is not seen in obser- 
vational data. The model also tends to under represent the 
abundance of massive (M, > 10 11 fr _1 M ) galaxies at z > 1 
compared to the B8W7 model and the observational data. 

• We also consider a model in which the wind mass- 
loading scales with the inverse of the circular velocity of the 
disk. This model results in a more subtle transition between 
the regime where material easily escapes the halo and that in 
which it stalls. With suitable choice of parameters, and the 
inclusion of the AGN feedback scheme, the model reproduces 
the observed galaxy mass function above M, ~ 10 9 /i -1 M Q . 
However, although the slope of the SSFR— M* relation is 
weaker than that seen in the SW model, the trend is re- 
versed compared to the slightly rising relation seen in obser- 
vational data. The model tends to under-predict the abun- 
dance of massive, z > 1 galaxies compared to B8W7. The 
deficit is, however, relatively small and the MS model could 
be compatible with the observational data if the random 
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uncertainties in the stellar mass estimates are greater than 
0.2 dex. 

• Feedback from AGN may have two very distinct forms. 
The "hot halo" (or "radio") mode feedback suppresses the 
supply of gas from cooling haloes. This is a key component 
of the successful models we present. However, it is interest- 
ing to investigate if this can be replaced by a "starburst" or 
"QSO" mode of feedback. We implement this by enhancing 
the wind speed and mass loading during starburst events. 
We find that while this mode may enhance the feedback from 
quiescent star formation, it does not introduce a character- 
istic break in the galaxy mass function and cannot be seen 
as an alternative to the "hot-halo" mode feedback. Similar 
conclusions are reached by Gabor et al. (2011), when they 
introduce a comparable feedback scheme in hydrodynamical 
models. 

• We have also investigated the star formation history 
of the each of the models, and its dependence on the stel- 
lar mass. Looking at the integrated star formation rate, the 
main difference lies in the strength of the decline between 
z — 2 and the present day. This is strongest in the B8W7 
model. Nevertheless, given the uncertainties in the obser- 
vational estimates, both the MS and B8W7 model seem to 
provide a reasonable description. However, the models show 
a much greater variety in behaviour when the results are 
broken down by stellar mass. For this reason, comparison 
of the total star formation rate may be miss-leading since 
the observational results assume large extrapolations to ac- 
count for galaxies that are not directly observed. It is there- 
fore cleaner to compare the contribution to the total star 
formation rate density at each epoch as a function of stel- 
lar mass. None of the models provide a perfect description 
of the data. At low redshift, the distribution is broader in 
the B8W7 model than is the case for the SW or MS mod- 
els. Observed contribution from low mass galaxies is better 
described by B8W7, but the model contains too few high 
star formation rate galaxies at low redshift. At z ~ 2, the 
situation is reversed, with all the models showing a similar 
contribution from low mass galaxies, and the primary differ- 
ence being the paucity of high stellar mass galaxies in the 
MS and SW models. 

• The differences between the models are emphasised by 
comparing the mass fractions of baryons in different phases 
as a function of halo mass. The conventional model, B8W7, 
assumes that haloes less massive than 10 ll s h~ 1 M.Q are 
baryonically closed. Almost all of the baryons in low mass 
haloes are placed in cool clouds in the halo (because of the 
short halo cooling time). In contrast the SW and MS, ex- 
pel much of their baryonic content from the halo. The frac- 
tion of baryons retained in the halo has a maximum of 80% 
at Mhaio = 1O 11,7 /i _1 M0 and drops rapidly to a minimum 
of ~ 10% at M ha i = lO n ft _1 M . At the minimum, the 
majority of baryons are predominantly stars and cold gas. 
At masses above 10 h~ Mq, the "hot-halo" mode of feed- 
back takes over expelling the halo material, resulting in low 
baryon fractions as discussed in Bow08, and in good agree- 
ment with X-ray observations of groups and clusters. The 
halo baryon fractions of the MS and SW models offer a good 
explanation of the low abundance of HI clouds around M31 
(Sancisi et al. 2008). Further work is required to compare 
the model predictions to such data and to the circumgalac- 
tic gas haloes inferred from absorption line systems at low 



redshift (Tumlinson et al. 2011) since the ionisation state of 
the extra-galactic gas must be carefully computed and com- 
bined with X-ray limits on the emission from galactic haloes 
(Crain et al. 2010). In all of these models, the recycling of 
previously ejected gas plays an important role. 

Another potential discriminant between the models is the 
distribution of metals in the ISM. All the models predict 
that a major fraction of the metals produced by stars will 
reside outside galaxies. In the absence of AGN feedback, 
the B8W7 predicts that the metals will be confined within 
galaxy haloes, while it will be more wide-spread in the ex- 
pulsion models (SW and MS). We find, however, that at 
z < 2 the major contributor of extra-galactic metals is the 
AGN powered hot-halo feedback that expels diffuse gas from 
galaxy groups. This material is highly enriched compared to 
the winds from low-mass galaxies: combined with the high 
fraction of the total baryon budget that has been expelled 
from groups, this component dominates the diffuse metal 
content of the low redshift Universe. At z > 3, AGN feed- 
back makes little contribution to the intergalactic fraction of 
metals and the observed wide-spread distribution of metals 
at this epoch favours the expulsion feedback schemes of the 
MS and SW models. However, the nature of our phenomeno- 
logical model makes it difficult to predict the dispersal of 
the expelled material, and it is unclear how far from the 
parent galaxy groups the metals will be spread (Booth et al. 
2010). Careful consideration of the dynamics of the outflow 
is needed to make a meaningful comparison with observa- 
tional data on metal absorption lines. 



In summary, although we have introduced a feedback 
scheme that reproduces the results of hydrodynamical sim- 
ulations well, we find that the original wind scheme of Bow06 
produces a better match to observational data on the stellar 
mass content of haloes, the specific star formation rates of 
galaxies and the evolution of these quantities. Comparison of 
the models with the observational data highlights several as- 
pects of the observational data that are not consistent with 
any of the models considered. Firstly, none of the models re- 
produce the rapid drop in the normalisation of the observed 
SMF between z = and z = 1 — 2. Secondly, although the 
B8W7 model reproduces the observed dependence of the 
specific star formation rate on stellar mass better than the 
SW or MS models, none reproduce the observed tendency for 
lower mass galaxies to have higher SSFR than their massive 
counter-parts. Both of these trends suggest that the effective 
wind speed (relative to the halo escape velocity) should in- 
crease with redshift. One possibility is that this might arise 
because of the more clumpy concentration of star forming 
regions in high redshift galaxies. This is an interesting possi- 
bility for further investigation. There is also no intrinsic rea- 
son why galaxy winds should be a simple power-law, and the 
observational suggest we should explore feedback schemes in 
which the exponent of the mass dependence varies with sys- 
tem. For example, the problem of the SW and MS models 
over producing galaxies below 10 9 F'Mg could be solved by 
introducing a strong mass dependence to the feedback below 
this limit. Rather than introducing phenomenological mod- 
ifications to the feedback scheme, an alternative approach 
would be to simulate idealised galaxies at high (10 pc), or 
ultra-high (O.lpc) resolution and to extract suitable param- 
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eterisations for the outflow Hopkins, Quataert, & Murray 
(2011); Stringer et al. (2011); Creasey et al. (2011). 

In future, we will combine the expulsion mode of feed- 
back considered here with developments of the GALFORM 
code (Lagos et al. 2010) to more accurately trace the cold 
gas content of the universe. We will also make a one-to- 
one comparison of the formation histories of galaxies from 
GIMIC, and the forth-coming suite of EAGLE simulations, 
with the GALFORM code. Our aim is not only to develop 
the GALFORM model as an emulator of numerical simu- 
lations (eg., Bower et al. (2010)) but also to use it as a 
fundamental tool for understanding the key components of 
successful galaxy formation models. Only by representing 
the evolution of galaxies as a simple set of coupled differen- 
tial equations can we claim to have understood the galaxy 
formation problem, and to have separated the key processes 
from the details. 
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